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PREFACE 


Work  on  a  program  to  study  the  effects  of  natural  phenomena  on  pre¬ 
cise  surveying  and  position  determination  methods  and  on  the  application  of 
natural  phenomena  to  geodesy  was  initiated  on  19  September  1958  under  Contract 
No.  DA-44-009-ENG-3769  with  the  Topographic  Systems  Research  Branch  of  the 
Topographic  Engineering  Department  of  the  U.  S.  Army  Engineer  and  Development 
Laboratories.  On  30  January  1960,  a  supplemental  agreement  was  executed  which 
extended  the  period  of  the  contract  to  19  December  1961  and  redirected  the 
activities  into  narrower  fields.  On  26  January  1962,  another  agreement  was 
executed  which  extended  the  period  of  the  contract  to  19  December  1963  and 
redirected  the  activities  to  the  study  of  a  particular  distance  measuring 
system.  The  work  conducted  after  19  December  1961  is  summarized  in  this 
report. 


During  the  current  contract,  the  activity  on  the  program  has  been 
under  the  direction  of  Mr.  P.  C.  Constant,  Head,  Electronics  Section;  Mr. 

R.  S.  Brown  has  been  the  project  leader,  and  is  responsible  for  the  theoret¬ 
ical  work.  Mr.  E.  J.  Martin,  Jr.,  has  been  responsible  for  the  experimental 
designs. 
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SUMMARY 


This  project  has  been  a  feasibility  study  to  determine  whether  it  was 
possible  to  deduce  the  index -of-refractlon  profile  along  some  path  in  the  atmos¬ 
phere  from  physical  measurements  made  external  to  that  path.  This  problem  was 
solved  and  it  was  demonstrated  that  it  was  possible  to  make  this  remote  meas¬ 
urement  of  the  index-of -refraction  profile  using  a  quantity  called  the  "tran¬ 
sient  reflection  coefficient"  of  the  atmosphere. 

When  an  electromagnetic  wave  propagates  into  a  medium  of  changing 
index-of -refraction,  a  certain  portion  of  the  Incident  energy  is  reflected  back 
toward  the  source.  The  ratio  of  the  reflected  energy  to  the  transmitted  energy 
is  called  the  reflection  coefficient  of  the  atmosphere.  It  was  first  shown  that 
this  reflection  coefficient,  which  is  a  function  of  the  frequency  of  the  inci¬ 
dent  wave,  is  uniquely  determined  by  the  index-of-refractlon  profile.  In  addi¬ 
tion,  it  was  demonstrated  that  if  the  value  of  the  reflection  coefficient  were 
known  for  all  frequencies,  then  it  would  be  possible  to  reconstruct  the  index - 
of-refraction  profile. 

Because  the  index-of-refraction  profile  changes  with  time,  the  reflec¬ 
tion  coefficient  is  also  a  time  variable.  Consequently,  in  order  to  determine 
the  value  of  the  reflection  coefficient  at  all  frequencies,  it  would  be  neces¬ 
sary  to  make  a  large  number  of  simultaneous  measurements  each  at  a  different 
frequency.  Such  a  system,  although  feasible,  would  not  be  practical.  It  was 
because  of  this  objection  that  the  concept  of  a  "transient  reflection  coeffi¬ 
cient"  was  introduced. 

When  a  very  short  pulse  is  transmitted  into  a  region  of  variable 
index  of  refractions,  the  reflected  energy  will  also  be  in  the  form  of  a  pulse 
rather  than  a  continuous  wave.  It  was  demonstrated  that  if  the  transmitted 
pulse  were  of  a  sufficiently  short  duration,  then  the  ratio  of  the  Fourier 
transforms  of  the  reflected  and  transmitted  pulses  would  contain  sufficient 
spectral  information  to  allow  the  construction  of  the  index-of -refraction  pro¬ 
file.  It  is  this  ratio  that  is  called  the  "transient  reflection  coefficient" 
which  differs  from  the  ordinary  reflection  coefficient  only  in  the  way  in  which 
it  is  measured.  The  principal  advantage  of  measuring  the  reflection  coefficient 
by  a  transient  method  is  that  only  a  single  measurement  need  be  made.  An  addi¬ 
tional  advantage  is  that,  because  the  ccoplete  measurement  can  be  made  in  a  very 
short  Interval  of  time,  the  adverse  effects  of  the  temporal  variations  of  the 
index-of-refraction  profile  are  nullified. 

Experiments  to  check  the  utility  of  the  inverse  scattering  theory  are 
proposed,  and  discussed.  These  experiments  consist  of  (a)  setting  up  a  trans¬ 
mission  path  whose  length  and  Index-of-refractlon  profile  are  accurately  known. 
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(b)  measurement  of  the  reflection  coefficient  of  the  transmission  path,  (c) 
calculation  of  the  index-of-refraction  profile  using  data  collected  in  (b), 

(d)  calculation  of  the  length  of  the  transmission  path  from  the  experimentally 
determined  index  of  refraction, and  (e)  the  results  of  (d)  will  be  compared  with 
a  true  length  of  the  line. 


I.  INTRODUCTION 


This  program  has  been  a  feasibility  study  of  a  new  technique  for 
using  the  transient  reflection  coefficient  of  the  troposphere  to  make  cor¬ 
rected  distance  measurements. 

Present  systems  for  making  precise  distance  measurements  transmit 
electromagnetic  energy,  such  as  light  or  radio  waves,  through  the  atmosphere 
to  a  reflector.  The  transmitted  energy  is  then  reflected  back  to  the  point 
of  transmission  and,  by  measuring  the  time  necessary  for  the  energy  to  make 
this  transit,  the  distance  between  the  transmitter  and  reflector  is  deduced. 

In  all  the  present  systems,  it  is  assumed  that  the  energy  travels  with  a  con¬ 
stant  velocity,  the  value  of  which  is  determined  by  the  local  atmospheric 
conditions  at  each  end  of  the  survey  path. 

The  index  of  refraction  varies  in  a  randan  manner  along  the  survey 
path.  These  variations  of  the  index  of  refraction  cause  corresponding  varia¬ 
tions  in  the  velocity  with  which  the  electromagnetic  energy  propagates,  and 
these  variations  cause  errors  in  the  indicated  path  length.  If  the  value  of 
the  index  of  refraction  were  known  everywhere  along  the  survey  path,  then 
these  errors  could  be  accounted  for  and  corrected.  The  task  of  this  program 
has  been  to  devise  a  method  for  measuring  this  index-of -refraction  profile. 

When  a  wave  propagates  into  a  medium  of  changing  index  of  refraction, 
not  only  does  the  velocity  of  the  wave  change,  but,  in  addition,  a  portion  of 
the  incident  energy  is  reflected  back  toward  the  source.  It  has  been  shown* 
that  the  index-of -refraction  profile  can  be  deduced  from  a  knowledge  of  the 
reflected  energy.  The  demonstration  of  this  fact  represented  the  successful 
completion  of  task  3a  as  specified  in  the  contract: 

"a.  Conduct  studies  and  investigations  to  determine  the 
feasibility  of  determining  the  index  of  refraction  of  the 
troposphere  between  two  points  to  be  used  for  a  correc¬ 

tion  to  a  distance  measurement  made  by  the  transmission  of 
electromagnetic  waves  between  the  two  points." 


*  Brown,  R.S.,  "Research  Studies  Related  to  Mapping,  Geodesy  and  Position 
Determination,"  14th  Interim  Technical  Report,  Contract  No.  DA -44-009- 
ENG-3769,  10  July  1962.  Also  see  Appendices  A  and  C,  this  report. 
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Subsequent  work  performed  under  this  contract  was  concerned  with 
task  3b,  the  second  and  final  provision  of  the  contract; 

"b.  Contingent  on  the  successful  accomplishment  of 
a.  above,  specify  the  system  characteristics  necessary  to 
determine  the  refractive  index  of  the  troposphere  between 
two  points,  e.g.,  transmitter  power,  pulse  width,  receiver 
sensitivity  and  noise  figure,  and  computational  circuitry." 

With  regard  to  task  3b,  only  a  qualitative  specification  of  the 
system  characteristics  could  be  made.  Any  further  efforts  toward  system 
specification  will  have  to  be  experimental  because  of  (a)  a  lack  of  knowledge 
of  the  characteristics  of  the  tropospheric  index-of -refraction  variations, 
and  (b)  the  necessity  of  extending  microwave  techniques  beyond  the  present 
state  of  the  art. 


II.  INVESTIGATION 


The  system  investigated  under  this  contract  is  shown  in  Fig.  la.  A 
radar  set  transmits  a  wave,  f(t)  ,  into  the  turbulent  region,  i.e.,  into  a 
region  of  variable  index  of  refraction.  As  the  transmitted  wave,  f(t)  ,  tra¬ 
verses  the  turbulent  region,  reflections  are  produced  and  this  reflected  wave, 
g(t)  ,  travels  back  to  the  radar  set  where  it  is  detected.  In  the  course  of 
this  investigation,  it  was  shewn  that  if  both  f(t)  and  g(t)  were  known, 
then  n(x)  ,  the  index-of -refraction  profile  characteristics  of  the  turbulent 
region,  could  be  uniquely  determined.  Now,  if  two  dielectric  markers  con¬ 
sisting  of  sheets  of  lucite,  teflon,  or  other  dielectric  material  were  intro¬ 
duced  into  the  propagation  path,  as  shown  in  Fig.  la,  then  they  would  produce 
two  characteristic  marks  in  the  index-of -refraction  profile  as  shewn  in  Fig. lb 
and  would  serve  to  mark  the  beginning  and  end  of  the  survey  path.  Thus,  the 
distance  between  the  dielectric  markers  would  be  determined  as  the  difference 
between  the  abscissal  values  of  the  dielectric  marks  shewn  on  the  index-of - 
refraction  profile. 


A.  Remote  Measurement  of  the  Index-of -Refraction  Profile 


A  wave,  w(k,x)  ,  propagating  in  a  medium  of  changing  index  of  re¬ 
fraction,  n(x)  ,  satisfies  the  differential  equation 


w(k,x)  +  l^n2(x)w(k,x)  •  0  ,  (l.l) 
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where  k  is  the  free  space  wave  number  of  the  wave  and  is  related  to  the 
temporal  frequency,  ia  ,  of  the  wave,  and  the  free  space  (in  vacuum)  velocity 
of  the  wave,  c  ,  by 


k  *  u)/c  .  (1.2) 

In  order  to  determine  n(x)  ,  the  equivalent  equation, 

u(k,x)  +  [k^-V(x)]  u(k,x)  »  0  (1.3) 

is  first  considered.  Here,  the  term  V(x)  ,  called  the  scattering  potential, 
is  equivalent  to  n(x)  . 

Now,  assume  that 


V(x)  »  0  for  x  <  0 


(1.4) 


and  that 


u(k,x)  =  eilDC  +  r(k)e"itoc  for  x  s  0  ,  (1.5) 


where  the  term  r(k)  ,  the  reflection  coefficient,  is  completely  known. 

In  Appendix  A  of  this  report,  it  is  shewn  that  r(k)  is  an  element 
of  a  certain  operator,  S  ,  called  the  "Scattering  Matrix": 


/■r(k)  T(k) 
^t(k)  p(k) 


(1.6) 


It  is  shown  that  the  elements  of  this  matrix,  which  are  functions  of  the  free 
space  wave  number  k  only,  are  uniquely  determined  by  the  scattering  potential, 
V(x)  . 
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Fortified  with  the  knowledge  of  the  uniqueness  of  S  ,  the  research 
efforts  were  then  directed  toward  the  solution  of  the  inverse  scattering  prob¬ 
lem,  namely,  the  problem  of  determining  V(x)  from  a  knowledge  of  S  .  This 
problem  was  solved  in  terms  of  a  function  K(x,y)  defined  by 


i_  j  „„  .  K(lt,y)  ♦  u.7) 

--  ItOOl2 


(Eq.  (A -41)  of  Appendix  A),  where  t(k)  is  the  transmission  coefficient 
(another  element  of  S  )  and  0(k, y )  is  a  certain  normalized  wave  (see  Eq. 
(A -38)  et  seq.  of  Appendix  A).  It  was  further  shewn  that  V(x)  and  K(x,y) 
were  related  by 


V(x)  •  2  J-  K(x,x)  (1.8) 

dx 


(see  Eq.  (A-55)  of  Appendix  A),  and  that  K(x,y)  could  be  found  by  solving 
the  integral  equation 

x 

0  »  R(x+y)  +  K(x,y)  +  £  R(y+z )K(x,Z )dz  ,  (1.9) 

-x 


where 


S(x.y)  .  i-  /’r(k).•11‘(x*J^,  dk 
2tt  » 


(1.10) 


(see  Eq.  (A-57)  et  seq.  of  Appendix  A). 

Equation  (1.9)  represents  a  solution  of  the  inverse  scattering  prob¬ 
lem.  It  was  shown  that  this  solution  would  be  valid  whenever  certain  condi¬ 
tions  pertaining  to  R(x+y)  were  met  (see  conditions  (l)  to  (5),  p.  42, 
Appendix  A  and  also  Eqs.  (A-60)  to  (A- 71),  Appendix  A). 
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The  next  task  was  to  determine  a  method  whereby  the  index-of -refrac¬ 
tion  profile,  n(x)  ,  could  be  constructed  from  the  scattering  potential  V(x)  . 
This  construction  is  given  in  Section  IV  of  Appendix  A  and  consists  of  the 
following  four  steps: 

1.  Find  U(z)  from 

U(z)  -  V(z )U(z )  =  0  ,  (1.11) 


2.  Find  x(z)  from 

x(z)  -  c  M5—  ,  (1.12) 

J  U2(z) 

3.  Invert  x(z)  to  obtain  z(x)  ,  and 

4.  Find  n(x)  from 

n(x)  ■  U2 (x )  .  (1.13) 


This  construction  completes  the  inverse  scattering  problem.  The  next 
task  was  to  investigate  methods  whereby  r(k)  could  be  measured. 


B.  The  Measurement  of  the  Reflection  Coefficient,  r(k) 

The  solution  of  the  inverse  scattering  problem  given  in  Section  II -A 
was  based  upon  the  assumption  that  r(k),  the  reflection  coefficient,  was 
known.  In  Appendix  B,  a  theory  is  presented  for  the  measurement  of  r(k)  , 
and  it  is  upon  this  theory  that  the  proposed  experimental  system  to  measure 
the  index  of  refraction  is  based. 

Very  briefly,  a  forward  traveling  wave,  F(x-ct)  ,  will  be  transmitted 
into  the  region  of  nonconstant  index  of  refraction  and  a  backward  traveling, 
reflected  wave,  0(x+ct)  ,  will  be  produced.  Both  of  these  waves  will  be  de¬ 
tected  and  recorded.  Then  the  quantities  f(k)  and  g(k)  will  be  computed. 
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(1.14) 


f  (k)  ■  -i—  f  F(x-ct)e"iIs(x"ct)d(x-ct) 

-NlSr  im 

and 

00 

g(k)  =  --=  f  G(x+ct)eik(x+ct)d(x+ct)  ,  (1.15) 

'v2TT  -ao 

and,  from  these,  r(k)  can  be  found  from 


r(k)  - 

f(k) 


(1.16) 


(see  Eqs.  (B-ll,  B-12,  B-13),  Appendix  B). 


C.  An  Alternate  Derivation  of  the  Inverse  Scattering  Theory 

The  representation  of  the  reflection  coefficient,  r(k)  ,  as  the 
ratio  of  the  spectrum  of  the  reflected  wave  to  that  of  the  transmitted  wave 
led  to  an  alternate  derivation  of  the  inverse  scattering  theory  in  which  the 
directly  measurable  terms,  f(k)  and  g(k)  ,  were  used  Instead  of  their  ratio. 
This  alternate  derivation  is  presented  in  Appendix  C,  and  follows  exactly  the 
same  course  as  the  first  derivation  except  that  the  wave,  u(k,x)  ,  was  ini¬ 
tially  normalized  as 


u(k,x)  *  f(k)eikx  +  g(k)e"il£X  for  x  s  0  ,  (1.17) 

Instead  of  as  in  Eq.  (1.5).  With  this  normalization,  K(x,y)  appears  at  the 
solution  of  the  integral  equation 

.x  x 

0  *  G(y+x)  +  J  K(x,z)F(y-z)dz  +  J  K(x,z)G(y+z)dz  (1.18) 

-X  -X 


instead  of 


0  -  R(x+y)  +  K(x,y)  +  J*  K(x,z )R(y+z )dz 

-x 
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(1.19) 


In  Eq.  (1.18),  the  functions  F(x)  and  G(x)  are  the  Fourier  transforms  of 
f(k)  and  g(k)  ,  respectively. 

The  two  inversion  formulae,  Eqs.  (1.16)  and  (1.19),  are  equivalent 
in  all  respects.  The  primary  difference  between  them  lies  in  the  fact  that 
Eq.  (1.18)  utilizes  directly  measurable  data,  whereas  the  function  R(x+y) 
appearing  in  Eq.  (1.19)  is  the  result  of  several  mathematical  operations  per¬ 
formed  on  F(x)  and  G(x)  . 

Two  examples  are  presented  in  Appendix  C  which  serve  to  indicate  the 
equivalence  of  the  two  inversion  formulae.  In  the  first  example,  the  func¬ 
tions 


and 


f(k)  -  iiisii) 

(kH)2 


G(k)  -  — - 

(k+i)3 


were  selected.  Then 


r(k) 


g(k)  „  -1 

f(k)  jj2+]i 


(1.20) 


(1.21) 


(1.22) 


It  was  then  shown  that  both  integral  equations  lead  to  the  same  differential 
equation, 


0  -  Kyy(x,y)  -  K(x,y)  +  K(x,-y)  ,  (1.23) 

whose  solution  is  (taking  into  account  the  boundary  conditions,  see  Eq.  (C-30), 
Appendix  C) 

K(x,y)  »  -  BlDhJ?x  +  .  (1.24) 

VsT  coshNSx 
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From  this  value  of  K(x,y)  ,  it  follows  that 


V(x)  =  -4  sinh2  */2x 


In  the  second  example,  the  functions 


and 


f(k) 


1 

k+i 


J 


g(k)  =-JL_ 

(k+i)2 


r(k)  = 

f(k) 


-i 

k+i 


were  selected.  In  this  case,  both  Integral  equations  lead  to  the 
equation 


0  *  Ky(x,y)  +  K(x,y)  +  K(x,-y)  , 


whose  only  solution  is 


K(x,y)  •  constant 


Thus,  in  the  second  example, 


V(x)  •  0  . 
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(1.25) 

(1.26) 

(1.27) 

(1.28) 

differential 

(1.29) 

(1.30) 

(1.31) 


D.  Calibrated  Transmission  Path 


To  check  the  utility  of  the  inverse  scattering  theory  Just  described, 
it  is  proposed  that  certain  physical  experiments  be  performed.  These  experi¬ 
ments  will  consist  of  setting  up  a  transmission  path  whose  length  and  index  - 
of -refraction  profile  are  accurately  known.  Next,  the  reflection  coefficient 
of  the  transmission  path  will  be  measured  and,  from  these  data,  the  index-of - 
refraction  profile  will  be  calculated  with  the  methods  of  either  Appendix  A 
or  Appendix  C.  Finally,  the  length  of  the  transmission  path  will  be  calculated 
from  the  experimentally  determined  index-of -refraction  profile,  and  this  re¬ 
sult  will  be  compared  with  the  true  length  of  the  line  to  obtain  an  estimate 
of  the  accuracy  which  can  be  expected  from  a  distance  measuring  system  based 
on  the  principles  presented  in  this  report. 

The  calibrated  transmission  path  will  consist  of  a  section  of 
band  waveguide  about  1  meter  long.  Associated  equipment  will  consist  of  a 
microwave  oscillator  covering  the  10  to  20  Gc.  band  of  frequencies  and  nec¬ 
essary  reflection  measuring  devices  (a  complete  listing  of  the  necessary  equip¬ 
ment  is  given  in  Appendix  D).  With  this  equipment,  the  length  of  the  line  can 
be  determined  to  an  accuracy  of  tl  mm.  (when  the  index  of  refraction  is  con¬ 
stant  ) . 


To  perturb  the  index-of -refraction  profile,  blocks  of  dielectric 
material  such  as  teflon  will  be  introduced  into  the  waveguide.  This  will 
cause  the  electrical  length  of  the  line  to  increase  and  will  also  cause  a 
certain  portion  of  the  electromagnetic  energy  propagating  within  the  line  to 
be  reflected.  These  reflections  can  be  measured  (on  a  steady-state  basis)  by 
standard  techniques,  far  instance,  with  a  slotted  line  or  with  a  Hewlett- 
Packard  4 ISA  ratio  meter. 

In  a  mare  practical  situation,  i.e.,  for  measuring  the  index-of - 
refraction  profile  of  a  turbulent  atmosphere,  these  steady-state  techniques 
will  fail  and  a  transient  method  must  be  used.  As  the  term  is  used  here,  a 
transient  measurement  is  one  in  which  the  complete  measure  of  the  index-of - 
refraction  profile  is  made  during  an  interval  which  is  short  with  respect  to 
the  rate  of  change  of  the  index- of -refract ion  profile.  To  do  this,  a  very 
short  impulse  of  electromagnetic  energy  will  be  transmitted  into  the  turbulent 
medium  and  a  corresponding  reflected  Impulse  will  be  received.  Then,  if  both 
the  transmitted  and  reflected  Impulses  are  known,  the  reflection  coefficient 
can  be  determined  with  the  methods  of  Appendix  B. 
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E.  Radio  Frequency  Energy  Source 


Methods  for  obtaining  the  short  pulses  (of  the  order  of  a  few  nano¬ 
seconds  in  duration)  of  radio  frequency  (rf )  energy  that  will  be  necessary  for 
the  actual  index-of -refraction  profile  measurements  have  been  investigated. 
Among  several  approaches  that  have  been  considered,  a  method  based  on  the  use 
of  ordinary  radar  pulse  techniques  with  subsequent  shaping  of  the  rf  pulse  has 
been  selected  as  most  promising.  This  method  contemplates  the  use  of  a  stand¬ 
ard  pulser,  a  magnetron,  and  a  rather  unconventional  application  of  a  gas- 
filled,  transmit-receiver  (T-R)  tube  as  a  pulse-shaping  device. 

Work  on  the  radio  frequency  energy  source  has  been  concerned  with 
the  specification  of  only  (a)  a  magnetron  and  a  suitable  pulser  (modulation 
equipment),  and  (b)  the  design  of  appropriate  control  and  monitoring  equipment. 
The  tentative  decision  to  perform  the  experiments  of  I^-band  frequencies  led 
to  the  selection  of  a  Type  QK  319  magnetron.  Past  experience  has  shown  that 
this  tube  can  be  successfully  driven  by  an  A PS-2  hard- tube  radar  modulator. 
Design  of  control  circuitry  has  included  means  for  selecting  any  one  of  nine 
crystal-controlled  pulse  repetition  frequencies  between  approximately  250  pps 
and  1,650  pps.  Means  for  selecting  any  me  of  three  pulse  durations  (0.5, 

1.0  and  2.0  p,sec.)  are  also  provided.  For  simplicity  and  flexibility,  it  is 
contemplated  that  the  radar  modulator  and  associated  control  circuitry  will 
be  mounted  in  a  single,  portable  unit.  This  portable  modulator  unit  will  also 
contain  a  Tektronix  RM  16  oscilloscope  for  monitoring  the  operation  of  the 
system. 


The  design  of  the  radio  frequency  energy  source  is  discussed  in 
detail  in  Appendix  E.  In  addition,  a  pulse  shaper  to  farm  nanosecond  pulses 
from  the  microsecond  pulses  generated  by  the  radio  frequency  energy  source 
will  be  required.  This  component  is  discussed  in  Section  III-B. 


III.  DISCUSSION 


The  problem  of  accurately  measuring  distances  is  of  fundamental 
importance  in  the  science  of  mapping,  geodesy,  and  position  determination. 
Variations  in  the  process  of  determining  the  position  of  a  point  "A"  with 
respect  to  another  point  "B"  can  be  obtained  through  application  of  different 
trigonometric  relations.  But  position  determination  must  ultimately  depend 
upon  either  distance  measurements  or  a  combination  of  distance  and  angle 
measurements . 


13 


For  centuries,  the  arts  of  surveying  and  mapping  were  completely 
dependent  upon  sene  form  of  the  level  or  transit  for  measuring  angles  and 
upon  seme  form  of  the  chain  or  tape  for  making  distance  measurements.  While 
proper  use  of  these  devices  can  result  in  highly  accurate  determination  of 
the  position  of  one  point  with  respect  to  another,  the  process  is  a  slow  one, 
especially  when  the  position  determination  is  to  be  carried  out  over  relatively 
large  distances  on  the  earth's  surface.  It  has  long  been  the  desire  of  sur¬ 
veyors  and  map  makers  to  be  able  to  accomplish  position  determination,  over 
long  distances,  rapidly  and  accurately.  The  development  of  radar  and  similar 
electronic  devices  during  World  War  II  provided  a  great  step  toward  the  reali¬ 
zation  of  this  goal. 

Since  1945,  a  number  of  military  equipments  have  been  tried  in  pre¬ 
cise,  long-distance,  measurement  experiments.  In  addition,  a  number  of  devices 
or  systems  have  been  developed  specifically  for  the  purpose  of  accurately 
measuring  long  distances.  Most  of  the  members  of  this  "first  generation"  of 
distance -measuring  instruments  depend  on  the  propagation  of  electromagnetic 
energy,  such  as  light  or  radio  waves,  through  the  atmosphere  to  a  reflector. 

A  portion  of  the  transmitted  energy  is  reflected  back  to  the  point  of  trans¬ 
mission,  and  the  distance  between  the  transmitter  and  the  reflector  is  deduced 
from  the  time  required  for  the  energy  to  make  this  transit. 

In  all  existing  systems  of  distance  measurement  by  means  of  electro¬ 
magnetic  propagation,  it  is  assumed  that  the  energy  travels  with  a  constant 
velocity,  the  value  of  which  is  determined  from  the  local  atmospheric  condi¬ 
tions  at  each  end  of  the  survey  path.  In  reality,  the  index  of  refraction 
can  vary  in  a  random  manner  along  the  survey  path.  Such  variations  of  index 
of  refraction  cause  corresponding  variations  in  the  velocity  of  propagation 
of  the  electromagnetic  wave.  These  variations  in  propagation  velocity  can 
produce  errors  in  the  indicated  path  length.  The  only  sure  means  of  elimina¬ 
ting  such  errors  depends  on  a  knowledge  of  the  index-of -refraction  profile 
along  the  entire  transmission  path. 

The  random  variations  of  the  index  of  refraction  cause  another 
effect,  namely,  the  reflection  of  a  part  of  the  Incident  electromagnetic 
energy  back  toward  the  source.  It  was  the  purpose  of  this  project  to  deter¬ 
mine  if  this  reflected  energy  could  be  utilized  to  obtain  more  precise  dis¬ 
tance  measurements.  The  major  result  of  the  project  was  the  theoretical 
demonstration  of  the  fact  that  the  reflected  energy  carried  with  it  sufficient 
information  to  construct  the  index-of -refraction  profile  along  the  survey  path. 
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A.  Determination  of  the  Index-of -Refraction  Profile 

The  starting  point  for  the  determination  of  the  index- of -refraction 
profile  was  the  Schrodinger  wave  equation 


u(k,x)  +  [k?-V(x)]  u(k,x)  =  0  .  (1.3)* 


In  this  equation,  the  term  u(k,x)  represents  a  wave  of  spatial  frequency  or 
wave  number  k  ,  propagating  in  the  x  direction.  The  term  V(x)  is  called 
the  scattering  potential  and  represents  the  variations  in  the  energy  (or 
velocity)  of  the  wave.  By  assuming  that  V(x)  *  0  for  x  <  0  ,  it  was  possi¬ 
ble  to  write  a  partial  solution  of  Eq.  (1.3 )  as 

u(k,x)  *  Aeilac  +  Be_ikx  ,  x  £  0  . 


In  this  representation,  the  term  A  exp(ikx)  represents  a  wave  propagating 
from  x  =  -«  into  the  region  where  V(x)  ^  0  and  the  term  B  exp(-ikx) 
represents  a  wave  reflected  fran  the  variations  of  V(x)  .  This  identifica¬ 
tion  is  consistent  with  the  assumption  of  a  suppressed  temporal  factor 
exp(-iuit)  . 


The  terms  A  and  B  represent  the  magnitude  of  the  transmitted  and 
reflected  waves,  respectively.  In  general,  both  A  and  B  are  functions  of 
frequency,  k  ,  and  are  related  by  the  equation 


|a|2  i  |b|2  , 


which  is  a  statement  of  the  lew  of  conservation  of  energy.  However,  because 
|A|2  is  a  measure  of  the  transmitted  power,  A  may  be  arbitrarily  selected 
as  being  A  ■  1  .  In  this  case,  the  partial  solution  of  the  Schrodinger 
equation  takes  the  special  form 

u(k,x)  ■  eilsx  +  r(k)e"ilot  ,  x  s  0  .  (1.5) 


*  Note  that  the  equation  numbers  used  in  this  section  are  identical  to  the 
numbers  used  on  the  corresponding  equations  in  Section  II. 
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In  Eq.  (1.5),  the  term  r(k)  Is  called  the  reflection  coefficient; 
it  is  a  measure  of  the  saount  of  energy  which  is  reflected.  It  is  also  an 
element  of  a  certain  matrix,  S  ,  called  the  scattering  matrix: 


In  Appendix  A,  the  significance  of  S  is  discussed  and  a  number  of  mathemati¬ 
cal  expressions  relating  its  elements  are  derived. 

Chce  the  scattering  matrix  and  its  relationship  to  the  scattering 
potential  had  been  defined,  it  was  possible  to  give  a  concise  mathematical 
statement  of  the  problem:  If  r(k)  is  known,  can  V(x)  be  found?  This 
problem  was  solved  by  presenting  a  method  whereby  V(x)  could  be  uniquely 
constructed  from  r(k)  . 

It  is  a  vell-kncwn  fact  that  both  u(k,x)  and  u(k,x)  are  con¬ 
tinuous  functions  of  x  .  Thus,  since  u(k,x)  was  known  for  all  negative 
values  of  x  ,  it  was  only  necessary  to  determine  a  method  whereby  this  known 
function  could  be  extrapolated  into  the  space  x  St  0  in  such  a  way  that 
u(k,x)  and  u(k,x)  would  be  continuous  and  such  that 


u(k,x)  -  [ V (x )-k^ Ju( k, x ) 


This  extrapolation  was  performed  by  Introducing  two  new  functions,  v(k,x) 
and  K(x,y)  .  The  first  of  these  functions,  v(k,x)  ,  was  defined  as  the  wave 
which  would  (a)  be  equal  to  u(k,x)  for  x  s  0  ,  and  (b)  would  satisfy  the 
Schrodinger  equation  for  x  at  0  if  V(x)  were  equal  to  zero.  The  second 
function,  K(x,y)  ,  was  not  explicitly  defined;  it  appeared  as  the  kernel  of 
the  integral  transform 


u(k,x)  -  v(k,x)  +  J  K(x,y)v(k,y)dy 
-x 


(a) 


(see  Eq.  (A -46),  Appendix  A).  The  use  of  this  extrapolation  formula  is  a 
rather  standard  technique  in  the  theory  of  partial  differential  equations. 

It  was  kneara  that  u(k,x)  obtained  from  Eq.  (a)  satisfied  all  of  the  required 
boundary  conditions  (this  fact  is  demonstrated  in  Appendix  A);  it  was  necessary, 


•  16 


however,  to  demonstrate  that  a  relationship  did  exist  (a)  between  K(x,y) 
and  V(x)  ,  and  (b)  between  K(x,y)  and  r(k)  .  This  demonstration  is  pre- 
sented  in  Appendix  A,  Eq.  (A -4 7 )  et  seq.  Very  briefly,  the  results  are : 

(a)  K(x,y)  is  related  to  V(x)  by 


V(x)  •  2  i-  K(x,x)  ,  (1.8) 

dx 


and  (b)  K(x,y)  is  related  to  r(k)  by 

x 

0  ®  R(x+y)  +  K(x,y)  +  J  R(y+z)K(x,z)dz  (1.9) 

-x 


where 


R(x+y)  -  i -  / r(k)e-lk(x+y)  dk  .  (1.10) 

2tt  v 


In  other  words,  if  r(k)  were  known,  then  K(x,y)  could  be  found  from 
Eq.  (1.9)  and  V(x)  from  Eq.  (1.6). 

Attention  is  now  turned  to  the  problem  of  determining  the  index-of- 
refraction  profile,  n(x)  .  If  u(k,x)  represents  an  electromagnetic  wave 
propagating  in  the  atmosphere,  then  it  satisfies  not  the  Schrodinger  equation, 
but  rather  the  Helmholtz  equation: 


u(k,x)  +  k?n2(x)u(k,x)  ■  0 


However,  the  effort  involved  in  finding  V(x)  was  worthwhile.  In  Appendix  A, 
Eq.  (A-72 )  et  seq.,  it  is  shown  that  if  V(x)  is  known,  then  the  correspond¬ 
ing  n(x)  can  be  found.  This  fact  is  demonstrated  by  a  standard,  but  rather 
involved,  change  of  variable  techniques;  the  reader  is  referred  to  Appendix  A, 
Eqs.(A-72)  to  (A-92)  for  the  details. 

The  next  task  in  this  study  was  concerned  with  the  method  whereby 
r(k)  could  be  measured.  This  problem  is  discussed  in  Appendix  B  where  it  is 
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shewn  that  if  a  short  impulse,  F(x-ct)  ,  is  propagated  into  the  medium  of 
varying  index  of  refraction  and  if  a  reflected  pulse,  0(x+ct)  ,  is  received, 
then  r(k)  can  be  found  freo  the  equation 


r 


(k)  =  fi i£l 
V  f(k) 


(1.16) 


where 


f (k)  *  — —  f  F(x-ct)e-ik^x"ct^d(x-ct)  ,  (1.14) 

•fin  -«> 


and 


g(k)  =  — —  /  G(x+ct)eik(x+ct  ^d(x+ct )  .  (1.15) 

•j  2n  _® 


B.  Proposed  Experimental  Verification  of  the  Theory 

On  the  basis  of  the  theoretical  results  obtained  under  this  program, 
it  is  possible  to  deduce  the  index-of -refraction  profile  from  measurements  of 
the  reflections.  However,  because  it  is  always  desirable  to  have  experimental 
evidence  to  support  a  theory,  an  experimental  research  program  to  obtain  these 
data  has  been  designed.  The  objectives  of  this  program  will  be: 

1.  To  verify  the  theory  presented  in  this  report; 

2.  To  design  and  assemble  an  experimental  system,  based  on  estab¬ 
lished  theory,  which  can  be  used  to  accurately  determine  the  index-of -refrac¬ 
tion  profile  along  a  prescribed  transmission  path;  and 

3.  To  develop  the  techniques  of  using  this  system  for  the  precise 
measurement  of  distance. 

To  achieve  these  objectives,  it  will  be  necessary  to  accomplish  the 
following  tasks,  in  essentially  the  order  listed: 
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1.  Establish  a  calibrated  transmission  path  for  electromagnetic 
waves,  along  which  the  index  of  refraction  can  be  precisely  controlled  and 
(within  limits)  changed  at  will. 

2.  Design  and  fabricate  a  radio  frequency  energy  source  capable  of 
generating,  on  a  repetitive  as  well  as  a  triggered  or  single-shot  basis,  very 
short  pulses  (of  the  order  of  a  few  nanoseconds  in  duration)  of  electromag¬ 
netic  energy  at  a  relatively  high  peak  power  level. 

3.  Design  and  fabricate  a  signal  separation  system  capable  of 
sampling  the  pulses  transmitted  by  the  system  described  in  (l)  and  the  re¬ 
flections  of  those  signals  from  the  calibrated  transmission  path  described 
in  (2). 


4.  Design  and  fabricate  a  monitor  system  capable  of  measuring  per¬ 
tinent  characteristics  of  the  direct  and  reflected  signals  produced  by  the 
signal  separation  system  described  in  (3). 

5.  Establish  data  processing  techniques  suitable  for  handling  the 
experimental  data  obtained  from  the  monitoring  system  described  in  (4), 
thereby  obtaining  the  index-of-refraction  profile  along  the  calibrated  trans¬ 
mission  path. 

6.  Perform  a  sufficient  number  of  experiments,  using  the  combina¬ 
tion  of  equipment  and  techniques  described  in  (l)  through  (5),  to  establish  the 
feasibility  of  accurately  measuring  distance  by  this  method. 

The  calibrated  transmission  path  mentioned  in  task  1  has  been  de¬ 
signed  and  the  major  components  are  listed  in  Appendix  D.  The  radio  frequency 
energy  source  (task  2 )  has  been  partially  designed  and  is  reported  in  Appendix 
E.  Further  design  must  be  made  on  an  experimental  or  cut-and-try  basis  because 
the  requirements  outlined  in  tasks  1  through  4  are  beyond  the  present  state- 
of-the-art  of  microwave  technology. 

The  proposed  research  program  could,  theoretically,  be  carried  out 
using  either  optical  or  radio  frequency  electromagnetic  energy.  However, 
there  are  certain  factors  which  make  microwaves  a  much  more  attractive  choice 
than  light  for  the  purposes  of  this  investigation.  First,  the  signal  that 
will  be  reflected  from  regions  of  varying  index  of  refraction  on  the  calibrated 
transmission  path  will  be  only  a  small  fraction  of  the  incident  signal;  pre¬ 
liminary  analysis  indicates  that  reflections  will  be  about  100  db  below  in¬ 
cident  signal  power  levels.  Consequently,  quantitative  reflected  signals 
would  be  extremely  difficult  to  measure  if  the  transmitted  signals  were  not 
incident  on  the  transmission  path  at  relatively  high  power  levels .  Recent 
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advances  in  laser  technology  have  resulted  in  devices  capable  of  producing 
well-collimated,  high-intensity  beams  of  light,  but  techniques  for  obtaining 
very  short  (nanosecond)  pulses  of  laser  output  have  not  been  fully  developed. 
Pending  advances  in  the  laser  field  may  alter  this  situation  in  the  near 
future,  but  microwave  techniques  which  offer  a  great  deal  of  premise  obtaining 
short-duration  pulses  of  rf  energy  are  in  existence.  Second,  even  if  short 
pulses  of  high-intensity  light  could  be  produced,  practical  limitations  on 
the  power  levels  of  incident  signals  that  could  be  handled  safely  in  the  pro¬ 
posed  experimental  system  would  result  in  reflected  signals  so  weak  as  to  be 
difficult  to  detect  by  any  known  optical  method.  This  is  not  believed  to  be 
the  case  for  micrcwave  signals  because  amplification  techniques  are  highly 
developed  in  the  field  of  microwave  electronics. 

On  the  basis  of  the  foregoing  considerations,  it  is  contemplated 
that  the  proposed  experimental  research  program  will  be  carried  out  using 
radio  frequency  (microwave)  electromagnetic  vaves.  The  exact  microwave  fre¬ 
quency  that  will  be  employed  depends  on  a  number  of  factors.  The  shorter 
wavelengths  of  the  higher  microwave  frequencies  offer  the  premise  of  greater 
inherent  accuracy  in  distance  measurement.  However,  the  variety  of  camer- 
cially  available  components  and  equipment  suitable  for  use  above  the  K-band 
frequencies  (18.0-26.5  Gc.)  is  very  limited.  Therefore,  K-band  is  considered 
optimum  for  the  proposed  program.  Other  frequencies  could  be  used,  with 
X-band  being  the  lower  limit. 

On  the  other  hand,  if  equipment  which  would  operate  in  the  millimeter 
wavelength  range  were  to  become  commercially  available,  it  would  be  interesting 
to  extend  the  range  of  investigation  up  to,  perhaps,  150  Gc.  to  evaluate  mare 
carefully  the  effects  of  atmospheric  attenuation.  A  curve  of  attenuation 
versus  frequency  is  shown  in  Fig.  2. 

Regardless  of  the  micrcwave  frequency  that  is  finally  selected  for 
the  performance  of  the  experiments,  the  general  nature  of  the  experimental 
system  will  be  as  shown  in  Fig.  3.  The  shaded  portions  of  this  figure  indi¬ 
cate  those  portions  of  the  system  that  can  be  easily  established  with  present- 
day  technology.  Although  seme  effort  will  be  required  to  procure,  assemble, 
and  adjust  these  portions  of  the  system,  no  major  technical  problems  are 
anticipated  in  these  areas.  The  major  part  of  the  research  effort  on  the 
proposed  program  will  be  devoted  to  the  design  and  assembly  of  those  portions 
of  the  system  shewn  by  the  unshaded  blocks:  (l)  the  mlse  shaper  in  the  rf 
energy  source,  (2)  the  signal  separation  system,  and  (3)  the  amplifier  in  the 
monitor  system. 

The  calibrated  transmission  path  will  consist  of  a  section  of  wave¬ 
guide  several  hundred  wavelengths  long  at  the  selected  operating  frequency. 
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ATTENUATION  (db/Km) 


Fig*  2  -  Attenuation  Characteristics  of  the  Atmosphere 
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Fig.  5  -  Block  Diagram  of  Proposed  Experimental  System 


An  additional  15-30  ft.  of  waveguide  will  be  used  to  produce  a  delay  in  the 
reflected  signal  so  that  switching  operations  can  be  performed  in  the  signal 
separation  system.  The  remainder  of  the  waveguide  path  will  be  sealed  and 
evacuated  to  provide  a  constant,  stable  atmosphere  in  which  the  actual  dis¬ 
tance  measurements  will  be  made.  However,  provisions  will  be  made  far  opening 
the  evacuated  line  so  that  dielectric  perturbations  of  known  characteristics 
can  be  introduced  into  the  waveguide  at  various  places.  The  entire  waveguide 
system  will  be  set  up  on  a  stable  foundation  which  will  include  provisions  for 
leveling  and  accurately  aligning  the  path. 

For  actual  index-of-refraction  profile  measurements,  the  tunable, 
c-w  sources  used  in  the  calibration  of  the  transmission  path  will  be  replaced 
by  a  pulsed  source.  This  pulsed  source  will  consist  of  a  magnetron  tube  and 
associated  modulation  circuitry  for  producing  pulses  of  rf  energy  either  on  a 
single-shot  or  repetitive  basis,  as  desired.  To  obtain  the  extremely  short 
pulses  (of  the  order  of  a  few  nanoseconds  in  duration)  necessary  for  the  pro¬ 
posed  experiments,  additional  microwave  circuitry  will  be  necessary  in  the 
form  of  a  pulse-shaping  system  at  the  output  of  the  magnetron.  The  develop¬ 
ment  of  this  pulse-shaping  system  represents  one  of  the  major  areas  of  re¬ 
search  effort  on  the  proposed  program. 

An  approach  to  the  solution  of  the  pulse-shaping  problem  is  suggested 
in  Fig.  4.  Here,  a  fast-ionizing  gas  tube  is  used  to  limit  the  transmitted 
signal  to  a  "leakage  spike"  which  1b  only  a  few  rf  cycles  in  duration;  a  fer¬ 
rite  circulator  device  provides  protection  for  the  magnetron.  The  use  of  a 
secondary-emission  "multipactor"  device  rather  than  a  gas  "transmit -receive" 
(T-R)  tube  presents  an  alternate  approach  to  the  solution  of  the  pulse-shaping 
problem.  Certain  new  types  of  voltage -controlled  ferrite  or  solid-state  wave¬ 
guide  switches  offer  additional  premising  approaches  to  the  solution  of  this 
problem. 


Since  the  determination  of  the  index-of-refractlon  profile  depends 
on  a  knowledge  of  the  signal  incident  upon  the  transmission  path,  as  well  as 
the  signal  reflected  from  variations  in  the  index  of  refraction  along  the 
path,  a  means  for  monitoring  both  incident  and  reflected  signals  is  an  essen¬ 
tial  portion  of  the  experimental  system.  The  fact  that  the  reflected  signal 
will,  very  likely,  be  quite  small  makes  the  separation  of  incident  and  reflec¬ 
ted  signals  rather  complicated.  Separation  of  the  signals  must  be  accomplished 
with  approximately  100  db  attenuation  of  the  incident  pulse  and  practically  no 
attenuation  of  the  reflected  signal.  The  development  of  such  a  signal  separa¬ 
tion  system  represents  a  second  major  area  of  research  on  the  proposed  program. 
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Fig.  4  -  Radio  Frequency  Pnlxe  Shaper 


A  possible  approach  to  the  solution  of  the  signal  separation  problem 
is  presented  in  Fig.  5.  The  ferrite  circulator  "1"  allows  the  signal  from  the 
pulse  shaper  to  be  transmitted  to  the  calibrated  transmission  path  with  prac¬ 
tically  no  attenuation  in  the  forward  direction  (the  direction  of  the  arrows). 
Leakage  in  the  reverse  direction  through  this  circulator  results  in  approxi¬ 
mately  20  db  of  attenuation  of  the  incident  signal  that  reaches  the  ferrite 
circulator  "2".  The  major  portion  of  the  signal  which  reaches  circulator  "2" 
is  transmitted  to  the  termination  "A".  Leakage  in  the  reverse  direction 
through  circulator  "2"  results  in  an  additional  20  db  of  attenuation  of  the 
incident  signal  which  reaches  circulator  "3".  The  process  is  repeated  through 
circulators  "4"  and  "5"  to  provide  approximately  100  db  of  attenuation  of  the 
incident  signal  which  is  ultimately  transmitted  to  the  monitor  system.  After 
a  sufficient  delay  to  allow  the  incident  pulse  to  reach  the  monitor  system,  a 
striking  voltage  is  applied  to  the  four  T-R  tubes  causing  them  to  ionize.  If 
the  section  of  waveguide  between  the  first  ferrite  circulator  and  the  cali¬ 
brated  transmission  path  is  sufficiently  long,  the  T-R  tubes  will  be  fully 
ionized  by  the  time  the  reflected  signal  returns  from  the  transmission  path, 
and  the  four  terminations  will  be  essentially  replaced  by  short  circuits.  The 
reflected  signal  will  be  readily  transmitted  from  the  circulator  "1"  to  the 
circulator  "2",  to  the  short  circuited  section,  back  to  circulator  "2",  to 
circulator  "3",  and  so  on  through  the  entire  signal  separation  system  (with 
only  a  few  decibels  of  attenuation)  to  the  monitor  system. 

Both  the  incident  and  reflected  signals  which  reach  the  monitor 
system  will  be  small.  Consequently,  the  first  portion  of  the  monitor  system 
must  provide  amplification  of  the  signal.  If  it  is  assumed  that  output  of  the 
pulse  shaper  comprises  1.0-nanosecand  pulse  of  50  kw.  peak  pcwer  at  1,000 
pulses/sec,  it  is  readily  seen  that  the  average  pcwer  output  of  the  rf  source 
is  about  50  mw.  Thus,  about  80  db  of  the  100  db  attenuation  introduced  by 
reflection  or  signal  separation  or  both  must  be  restored  by  an  appropriate 
amplifier  if  a  conveniently  detectable  0.5  aw.  signal  level  is  to  be  maintained 
at  the  detector  system  input.  The  assembly  of  a  suitable  amplifying  system  is 
the  third  major  area  of  research  effort  associated  with  the  proposed  program. 
Traveling  wave  tubes  and  parametric  amplifiers  offer  premise  in  meeting  the 
requirements  for  restoration  of  the  100  db  attenuation. 

The  remainder  of  the  monitor  system  comprises  instrumentation  for 
measuring  pertinent  characteristics  of  both  the  transmitted  and  reflected 
signals .  Amplitude  versus  time  information,  obtained  with  a  nanosecond  sam¬ 
pling  oscilloscope,  and  power  level  measurements,  obtained  with  a  bolometer 
and  rf  pcwer  meter,  will  be  necessary  for  the  determination  of  the  index-of- 
refraction  profile  along  the  path.  It  is  contemplated  that  coanerclally 
available  instruments  will  suffice  for  the  collection  of  these  data,  at  least 
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in  early  phases  of  the  program.  Modifications  of  conmercial  equipment  or  the 
development  of  specialized  instrumentation  may  become  necessary  in  later  phases 
of  the  program. 

Processing  of  experimental  data  will  be  accomplished  manually  or 
with  the  help  of  Midwest  Research  Institute's  computer  facilities.  Although 
automatic  data  processing  will  be  desired  at  some  future  time,  practically  no 
effort  in  the  development  of  a  special  data  processing  system  is  contemplated 
in  connection  with  this  program.  Since  the  data  processing  requirements  are 
well  defined  by  the  theory  which  forms  the  basis  of  the  proposed  experimental 
program,  only  a  small  amount  of  project  effort  will  be  necessary  to  establish 
the  computational  procedure. 

Experiments  will  actually  begin  with  the  testing  of  the  various 
systems  and  combinations  of  systems  as  they  are  developed.  However,  once  the 
over-all  experimental  system  has  been  perfected,  efforts  will  be  directed 
toward  the  perfection  of  techniques  for  operation  of  the  system  in  the  precise 
measurement  of  distance.  Also, experimental  evidence  of  the  system's  over-all 
capabilities  will  be  gathered. 


IV.  CONCLUSIONS 


1.  The  reflection  coefficient,  r(k)  ,  associated  with  an  index-of- 
refraction  profile,  n(x)  ,  can  be  determined  from  the  Bpectral  interpretation 
of  two  impulses  of  electromagnetic  energy. 

2.  If  r(k)  is  known,  then  n(x)  can  be  computed. 

3.  Experimental  verification  of  the  theory  presented  in  this  report 
should  be  pursued. 


V.  OTHER  APPLICATIONS  OF  THE  THEORY 


Although  the  theory  presented  in  this  report  was  developed  for  the 
purpose  of  obtaining  highly  accurate  distance  measurements,  it  has  a  number  of 
other  far-reaching  applications.  Basically  the  theory  indicates  the  method 
whereby  the  characteristics  of  a  medium  in  which  a  wave  can  propagate  can  be 
measured  remotely  and  rapidly  and  without  appreciably  disturbing  the  medium. 
For  instance,  not  only  could  the  spatial  characteristics  of  the  index  of 
refraction  be  measured  but,  by  making  a  sequence  of  such  measurements,  the 
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temporal  structure  of  the  turbulence  could  also  be  investigated.  An  investi¬ 
gation  such  as  this  would  provide  data  invaluable  to,  for  instance,  the  study 
of  energy  transfer  in  turbulent  flow. 

Other  applications  of  the  theory  arise  in  the  study  of  plasmas.  If 
a  transverse  electric  wave  penetrates  a  plasma,  then  its  behavior  within  the 
plasma  is  described  by  the  equation 


u(k,x) 


4ne^qj^N(x) 

m(cu^-v2) 


u(k,x)  =  0 


(5.1) 


where  N(x)  represents  the  electron  density  of  the  medium,  v  the  collision 
frequency  of  the  electrons,  <u  the  frequency  of  the  incident  wave,  e  and  m 
the  charge  and  mass  of  an  electron,  respectively,  and  all  other  symbols  are 
as  previously  defined.  For  very  high  frequencies  or  in  the  case  of  a  tenuous 
plasma,  v2  is  very  much  less  than  u£  and  may  be  neglected.  Thus,  Eq.  (5.1) 
becomes 

u(k,x)  +  ^k2  -  N(x)J  u(k,x)  *  0  .  (5.2) 


Equation  (5.2)  is  exactly  the  same  as  the  Schrodinger  equation  considered 
previously  with  the  term  V(x)  replaced  by 


V(x)  *  n(x) 
m 


(5.3) 


Thus,  the  electron  density  of  a  plasma  may  be  deduced  from  the  characteristics 
of  the  reflections  produced,  and  it  may  be  concluded  that  the  theory  presented 
here  is  applicable  to  the  investigation  of  plasmas,  including  shock  waves  and 
re-entry  sheaths. 

Waves  propagating  in  a  magnetoionic  medium  also  obey  Eq.  (5.2) 
whenever  the  electric  vector  of  the  wave  lies  parallel  to  the  direction  of  the 
superimposed  magnetic  field.  For  this  reason,  the  theory  is  also  applicable 
to  the  design  of  an  ionosphere  probe.  This  would  be  a  particularly  important 
application  because,  using  this  theory,  it  would  be  possible  to  make  a  measure¬ 
ment  in  an  interval  of  time  that  was  short  compared  to  the  rate  of  change  of 
the  ionospheric  structure;  consequently,  information  concerned  with  the  space- 
time  development  of  the  ionosphere  could  be  Obtained. 


28 


APPENDIX  A 


REMOTE  MEASUREMENT  OF  THE  INDEX-OF -REDACTION  PROFILE 
MATHEMATICAL  THEORY 


I.  Introduction 


The  problem  to  be  considered  here  is  that  of  determining  the  index- 
of -refraction  profile  along  some  prescribed  path  from  measurements  made  ex¬ 
ternal  to  that  path.  To  make  this  determination,  a  sharp  impulse  of  electro¬ 
magnetic  energy  will  be  propagated  into  a  medium  whose  index  of  refraction 
changes  from  point  to  point.  As  the  impulse  traverses  the  medium,  the  energy 
stored  in  its  electric  and  magnetic  fields  will  be  continuously  redistributed, 
and  a  portion  of  the  energy  will  be  reflected  back  toward  the  source.  In  the 
following,  it  will  be  shown  first,  that  the  characteristics  of  the  reflected 
impulse  are  uniquely  determined  by  the  index-of -refraction  profile,  and 
second,  that  the  index-of -refraction  profile  can  be  determined  from  measure¬ 
ments  of  the  reflection.  Prom  a  physical  point  of  view,  this  fact  represents 
a  possible  solution  to  the  problem  because  the  reflected  impulse  can  be  meas¬ 
ured  at  an  arbitrary  point  located  external  to  the  path. 

Any  wave,  w(k,x)  ,  propagating  in  a  medium  with  index  of  refraction 
u(x)  satisfies  the  Maxwell  wave  equation 


w(k,x)  +  k?n2(x)w(k,x)  *  0 


where  k  is  the  wave  number  or  spatial  frequency  of  the  wave.  Far  reasons 
of  mathematical  expediency,  however,  we  shall  at  first  concern  ourselves  not 
with  the  Maxwell  wave  equation,  but  with  the  Schrodinger  wave  equation 


u(k,x)  +  ^-Vfx)]  u(k,x) 


0 


where  V(x),  called  the  scattering  potential,  is  analogous  to  the  index- of - 
refraction  profile.  These  two  equations  can  be  tied  together  by  assuming  that 
n(x)  *  1  and  that  V(x)  «  0  for  x  <  0  .  Thus,  the  two  equations  are  identi¬ 
cal  for  x  <  0  and  we  can  assume  that  u()^x)*  w(k,x)  far  x  <  0  . 
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Subsequently,  In  Section  II  of  this  Appendix,  It  will  be  shewn  that 
the  Incident  and  reflected  waves  are  related  by  a  certain  operator  S  ,  called 
the  scattering  matrix.  Study  of  this  operator  will  yield  considerable  Insight 
to  the  nature  of  the  reflection  coefficient.  Then,  in  Section  III,  it  will  be 
shown  that  the  scattering  potential,  V(x)  ,  can  be  determined  if  the  reflec¬ 
tion  coefficient  is  known.  In  Section  IV  it  will  be  shown  how  n(x)  can  be 
determined  from  V(x)  . 


II.  The  Scattering  Matrix 

Scattering,  or  radiative  transfer,  processes  are  of  fundamental 
importance  in  many  fields  of  mathematical  physics  such  as  electromagnetic 
theory,  quantum  mechanics,  neutron  diffusion,  etc.  Usually,  the  aim  is  to 
determine  the  asymptotic  scattered  field  in  terms  of  the  incident  field  and 
some  prescribed  potential  field  or  medium  which  is  responsible  for  the 
scattering.  However,  in  this  study,  interest  will  be  focused  on  the  inverse 
problem,  that  of  specifying  the  characteristics  of  the  scatterer  in  terms  of 
the  incident  and  scattered  fields. 

This  discussion  will  be  concerned  with  the  one-dimensional  case 
which  is  governed  by  the  Scbrodinger  wave  equation 


ii(k,x)  +  £k^-V(x)  j  u(k,x)  =  0 


(A-l) 


and  it  will  be  assumed  that 


V(x)  *  0  for  x  <  0 


(A-2) 


throughout  this  discussion,  although  the  results  to  be  obtained  will  be  valid, 
or  can  be  easily  generalized,  when  a  less  severe  condition  replaces  Eq.  (A-2). 
For  purposes  of  illustration,  it  will  also  be  assumed  that 


V(x)  *  0  f or  x  >  x0  .  (A-3) 

However,  it  will  become  apparent  (because  xQ  nowhere  enters  into  the  argu¬ 
ment)  that  this  condition  is  invoked  purely  for  purposes  of  illustration  and 
represents  no  limitation  to  the  theory. 
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Let  u(k,x)  be  any  solution  of  Eq.  (A-l).  Then  it  can  be  decomposed 
into  two  parts,  the  incident  and  the  reflected  wave: 

u(k,x)  *  Ui(k,x)  +  ur(k,x)  .  (A-4) 

Consider  first  the  incident  wave  u^(k,x)  ;  it  can  be  further  decomposed  into 
two  parts,  one  incident  on  the  perturbation  from  the  left: 

u^kjx)  =  A_eikx  for  x<0  (A-5) 


and  one  incident  from  the  right: 


u^(k, x)  =  A+e'ikx 


for  x  >  xQ 


(A-6) 


Let  it  be  emphasized  at  this  point  that  the  only  significance  of  the  condition 
implied  by  Eq.  (A-3)  is  that  it  allows  one  to  speak  of  the  "right  side  of  the 
perturbation";  if  that  condition  had  not  been  invoked,  it  would  be  necessary 
to  write  Eq.  (A-6)  as 


lim  [u(k,x)-A+e’itat]  *  0  .  (A-7) 

x-*» 


In  a  similar  manner,  the  reflected  wave  can  be  decomposed  into  two  components, 
a  right-hand  component  and  a  left-hand  component : 


••iloc 

Uj.(k,x)  *  B_e  for  x  <  0  (A-8) 

and 

iloc 

Up(k,x)  ■  B+e  for  x  >  x0  .  (A-9) 


31 


Combining  the  above  notation. 


u(k,x)  =  A_eikx  +  B.e"lkx  for  x  <  0  (A-10) 

and 

u(k,x)  =  A+e'iloc  +  B+eiloc  for  x  >  xQ  .  (A-ll) 

It  is  convenient  to  consider  the  incident  wave  components,  A.  and 
A+  ,  as  being  "prescribed"  quantities  and  the  reflected  or  scattered  compo¬ 
nents,  B_  and  B+  ,  as  being  "measured"  quantities.  Clearly,  B+  depends 
in  a  linear  fashion  upon  both  A+  and  upon  V(x)  .  This  dependence  can  be 
easily  expressed  in  terms  of  an  operator  S  ,  the  scattering  matrix,  a  nota¬ 
tion  introduced  by  Heisenberg  in  1943: 


( 


B. 

B+ 


) 


(A-12) 


This  matrix  S  is  a  2  x  2  square  matrix  having  the  elements 


/r(k)  r(k)\ 

S  ■(  )  .  (A-13 ) 

\t(k)  p(k)/ 


Clearly,  these  elements  are  determined  by  V(x)  alone  and  are  functions  of 
the  complex  variable  k  .  Before  proceeding  to  the  actual  discussion  of  the 
inverse  problem,  it  is  expedient  to  derive  sane  relations  concerning  the 
elements  of  S  . 

The  first  item  of  interest  is  the  fact  that 

kl2  +  kl2  -  |b_|2  ♦  k|2  .  (A-i4 ) 
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Equation  (A-14)  is  no  more  than  a  statement  of  the  first  lav  of  thermodynamics. 
The  term  on  the  left  of  Eq.  (A-14)  is  a  measure  of  the  energy  incident  on  the 
scatterer;  the  term  on  the  right  is  a  measure  of  the  scattered  energy;  and 
because  V(x)  is  real,  the  scattering  process  is  lossless  and  the  equality 
must  hold.  From  a  mathematical  point  of  view,  Eq.  (A-14)  can  be  proved  quite 
easily  by  considering  the  Wronskian  of  tvo  solutions  of  Eq.  (A-l). 

If  u  and  v  are  tvo  solutions  of  Eq.  (A-l)  (or,  for  that  matter, 
of  any  differential  equation),  then  W(u,v)  ,  the  Wronskian  of  these  solutions, 
is  defined  as  the  determinant 


W(u,v) 


u  v 

u  v 


(A-1S ) 


vhere  the  dots  indicate  differentiation  vith  respect  to  x  .  It  vill  first 
be  demonstrated  that  W(u,v)  is  a  constant.  Differentiate  W(u,v)  vith 
respect  to  x  to  obtain 


—  W(u,v) 
dx 


U  V 

•  •  #• 
U  V 


(A-16) 


Substitute  the  values  of  u  and  v  into  Eq.  (A-16): 


—  W(u,v) 
dx 


1  1 
[k?-V(x)]  [k^-V(x)] 


(A-17) 


Then 

££  W(u,v)  *  0 
and 


W(u,v)  *  constant 


(A -18 ) 
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Now  to  prove  Eq.  (A-14)  let  u  be  the  solution  given  by  Eqs.  (A-10)  and 
(A-ll)  and  let  v  be  the  complex  conjugate  of  u  .  Then, 


and 


W  »  21k 


B_  A* 
A.  B* 


for  x  <  0 


W  -  2ik 


a+  b; 

b+  a: 


for  x  >  xQ 


(A-19) 


(A-20) 


Because  W  is  constant,  these  two  expressions  can  be  equated  to  obtain 
Eq.  (A-14). 

It  will  next  be  shown  that  S  is  unitary;  this  follows  immediately 

from  the  fact  that  (a)  S  transforms  the  vector  (^”)  into  the  vector  ^ 

(see  Eq.  (A-12)),  (b)  the  length  of  these  two  vectors  is  equal  (see  Eq.  (A-14)), 
and  (c)  any  operator  which  transforms  a  vector,  having  a  finite  number  of 
components,  into  another  vector  of  equal  length  is  unitary.  An  algebraic 
proof  of  this  statement  can  be  obtained  by  substituting  values  of  B_  and 
B+  ,  obtained  from  Eq.  (A-12),  into  Eq.  (A-14)  to  obtain 


1A-I2[1-Ir|2-It|2]+  |a+|2[i-M2-|p|2]-  A_A?[rT*+tp*] 

-  A2A+[r*T+t#p]  *  0  .  (A-21) 

Because  A_  and  A+  are  arbitrary,  it  is  necessary  that 


(A-22) 

(A-23) 

rr *  +  tp*  *  0  . 

(A “24 ) 
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By  making  use  of  these  last  three  expressions,  It  follows  immediately  that 


S' 


1 


(A-25 ) 


which  proves  the  assertion. 

The  last  result  of  interest  at  this  time  is  the  reciprocity  law: 


t  »  T 


(A-26 ) 


To  prove  the  reciprocity  relation,  first  write  down  the  matrix  representation 
of  u*(k,x )  remembering  that,  in  the  conjugate  wave,  the  roles  of  the  inci¬ 
dent  and  scattered  waves  are  interchanged: 


A*  «  SB* 


(A-27) 


Then 


B*  •=  S*ft* 


(A-28) 


However,  the  complex  conjugate  of  Eq.  (A-12)  is 


B*  «  S*A*  .  (A-29 ) 

Comparison  of  Eq.  (A-28)  with  Eq.  (A-29)  shows  that 


S  *  Sip 


(A -60) 


which  proves  the  assertion. 
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As  mentioned  previously,  the  vector  can  be  considered  to  be 

arbitrary.  Thus,  by  selecting  A_  =  0  and  A+  «  1  ,  ve  obtain 

u(k,x)  ■  eikx  +  r(k)e‘itoc  ,  for  x  <  0  (A-31) 

and 

lim  [u(k,x)-t(k)eikx]  =  0  .  (A-32) 

X— 

III-  The  Inverse  Scattering  Problem 

The  wave  u(k,x)  is  assumed  to  satisfy  the  Schrodlnger  wave  equation 

vi(k,x)  +  [k2-V(x)]  u(k,x)  *  0  (A-33) 

where  the  scattering  potential,  V(x)  ,  is  a  real,  piecewise  continuous  function 
of  x  and 

V(x)  =  0  for  x  <  0  .  (A -34 ) 

Consequently,  u(k,x)  can  be  defined  for  negative  values  of  x  as 

u(k,x)  «■  eikx  +  r(k)e'ikx  ,  x  <  0  .  (A-35) 

The  inverse  scattering  problem  is  that  of  constructing  an  expression  for  V(x), 
valid  for  all  x  ,  frcm  the  known  reflection  coefficient,  r(k)  . 

let  v(k,x)  be  a  wave  solution  of  the  differential  equation 

v(k,x)  +  k^C^x)  -  0  (A -36 ) 
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vhich  is  equal  to  u(k,x)  for  x  <  0  .  By  the  term  "wave  solution",  we  mean 
a  solution  whose  first  derivative  exists  and  is  continuous  for  all  values  of 
x  .  Consequently, 


v(k,x)  *  eilsx  +  r(k)e'iloc  ,  all  x 


(A-37) 


Let 


v(k,x) 


be  another  wave  solution  of  Eq.  (A-36)  and  be  given  by 


v(k,x)  =  e'1*®  -  r*(k)eiloc 


(A -38) 


where,  as  usual,  the  asterisk  indicates  a  complex  conjugate  quantity.  These 
two  wave  solutions  were  selected  because  they  have  the  property  that 

f  vfcxftfay.)  ^  =  6(x.y)  (A-39 ) 

2x1  -co  |t(k)|2 

where  t(k)  is  the  transmission  coefficient  associated  with  r(k)  and 
6(x-y)  is  the  Dirac  delta  function. 


»<x-y)  -  .“<*•*>  dk 

-OB 


(A-40) 


From  another  point  of  view,  the  function 


|t(k)l2 


may  be  looked  upon  as  the 


kernel  of  a  certain  integral  transform  acting  on  v(k,x)  .  Noting  the  simi¬ 
larity  between  Eqs.  (A -33)  and  (A-36),  we  apply  the  same  integral  transform 
to  u(k,x)  and  write  the  result  in  terms  of  a  new  function  K(x,y): 


dk  -  K(x,y)  +  6(x-y)  .  (A-41) 

It(k)|2 

Now,  if  this  new  function  K(x,y)  does,  in  fact,  exist  (we  shall  prove 
its  existence  later),  it  is  easy  to  shew  that 
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(A-42) 


±-  f  f  djdk  -  v(k,x)  +  f  K(x,z)v(k,z)dz 

2n  V-a  J-w  ! t ( J  )|2  ‘too 


whereas 


—  f  f  sLAji^Slaidial  ^  .  u(k>x)  .  (A^3) 

*"L  |t(j)l2 

New,  under  very  general  conditions  (Fubini's  theorem)  Eqs.  (A-42)  and  (A -43) 
can  be  equated  to  give 

00 

u(k,x)  =  v(k,x)  +  /  KfozMkjzJdz  .  (A-44) 


However,  because  u(k,x)  *  v(k,x)  for  x  <  0  ,  K(x,y)  is  triangular  in  the 
sense  that 

K(x,y)  *  0  far  y  <  -x  .  (A-45) 

Thus,  Eq.  (A-44)  becomes 

x 

u(k,x)  *  v(k,x)  +  K(x,z)v(k,z)dz  .  (A-46) 


It  is  now  necessary  to  prove  the  existence  of  K(x,y)  .  To  do  this, 


apply  the  operator 


d2  1 

-  -  V(x)  to  Eq.  (A-46).  After  substitution  from 

.dx2  J 


Eqs.  (A -33)  and  (A-36),  this  gives 


38  - 


(A-47 ) 


-k2u(k,x)  =  -k2v(k,x)  -  V(x)v(k,x)  -  V(x)  J  K(x,y)v(k, y)dy 


d2 
dx2  l 


X 

f  K(x,y)v(k,y)dy 


The  derivative  in  Eq.  (A-47)  is  evaluated  according  to  Leibnitz'  rule  to  yield 


-k2u(k,x)  =  -k2v(k,x)  -  V(x)v(k,x)  +  v(k,x) 

dx 


+  K(x,x)  — +  v(k,-x)  +  K(x,  -x)  ^  v(k, -x) 


Kx(x.,x)v(k,x)  +  l^(x,-x)v(k,-x)  +  P  K^foyM^yJdy 

-X 


A 

-  V(x)  J  K(x,y)v(k,y)dy 


(A -48 ) 


Next,  Eq.  (A-46)  is  multiplied  by  -k2 


A 

-k?u(k,x)  *  -J^vfk, x )  +  /  K(x,y)  [-l^vfoy)]  dy  .  (A-49) 

-x 


The  bracketed  term  in  the  integral  of  Eq.  (A-49)  is  evaluated  from  Eq.  (A-36). 
Thus,  frcm  Eq.  (A-49) 


x  2 

-k2u(k,x)  ■  -k2v(k,x)  +  P  K(x,y)  ~  v(k,y)dy  .  (A-50) 

•/„  dy2 
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The  integration  indicated  in  Eq.  (A-50)  can  be  performed  twice  by  parts  to 
yield 

-k2u(k,x)  »  -k2v(k,x)  +  K(x,x)  —  ^ 

+  K(x,-x)  -  Ky(x,x)v(k,x) 

+  Ky(x,-x)v(k,-x)  +  J  Kyy(x,y)v(k,y)dy 

-X 

(A-51) 

Next,  Eq.  (A-5l)  is  subtracted  from  Eq.  (A-48)  and  use  is 

relation 

made  of  the 

K  (x,x)  ♦  K  (x,x)  = 

*  y  dx 

(A-52) 

0  “  /  |^at(x,y)-i^y(x,y)-v(x)K(x,y)|v(k,y)dy 

+  ^2  -V(x)J  v(k,x)  +  2  dK^~X^  v(k,-x) 

(A-53 ) 

Equation  (A-53)  will  be  satisfied  if  the  following  three  conditions 

are  valid: 

KXx(X'y>  '  Kyy(X,y)  "  V(x)K(x'y)  "  0  ' 

(A -54) 

2  tf&LSi  -  V(x)  , 
dx 

(A-55) 

and 

axfrj-x)  . 
dx 

(A-56  ) 
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Now,  we  know  from  the  theory  of  second  order  differential  equations  that  a 
solution  of  Eq.  (A-54)  which  satisfies  the  boundary  conditions,  Eqs.  (A-55) 
and  (A -56 ),  does  exist.  Thus,  ve  have  proved  the  existence  of  K(x,y)  and 
nay  assume  Eq.  (A-46)  to  be  valid. 

The  next  step  is  to  find  K(x,y)  .  To  this  end,  apply  the  operator 


dy 


(A-57)* 


to  Eq.  (A-46)  (c  is  a  contour,  extending  from  K  *  -®  to  K  *  and  lying 
in  the  upper  half  of  the  complex  K- plane  above  any  poles  of  r(k)  ).  This 
yields 


0  =  R(x+y)  +  K(x,y) 


A 

f  R(y+z)K(x,z)dz 
J-x 


(A-58) 


where 


R(x+y)  =  —  fr(  k)e"ik^x+y^dk  . 
2t r  J 


(A-59) 


Equation  (A-58)  is  a  Fredholm  integral  equation  of  the  second  kind  and  can  be 
solved  for  K(x,y)  .  Thus,  Eqs.  (A-55),  (A-58),  and  (A-59)  represent  a  solu¬ 
tion  to  the  inverse  scattering  problem. 

The  above  presentation  represents  only  a  theoretical  solution  of 
the  inverse  scattering  problem.  From  a  practical  point  of  view,  the  solution 
is  not  yet  complete.  In  order  to  apply  the  theory  so  far  expounded,  it  is 
necessary  that  r(k)  be  known  for  every  value  of  k  .  However,  because 
r(k)  vill  be  determined  experimentally,  it  can  at  best  be  known  on  a  bounded 
real  interval  or,  even  more  likely,  for  only  a  finite  set  of  discrete  fre¬ 
quencies.  Hew,  because  not  every  function  r(k)  can  be  used  to  construct  a 
potential  v(x)  ,  a  very  real  problem  does  exist.  Accordingly,  we  wish  to 
determine  those  requirements  which,  when  placed  on  r(k)  ,  are  sufficient  to 
insure  that  a  scattering  potential,  v(x)  ,  can  be  constructed. 


*  For  instance,  this  operator  applied  to  the  function  f(y)  would  be 
written  i-  J  f(y)e“i^rdy  . 
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First,  let  us  consider  R(x)  .  Ve  state,  and  shall  subsequently 
prove,  that  the  following  requirements  are  sufficient  to  insure  that  V(x) 
can  be  constructed  from  Eq.  (A-55): 

1.  R(x)  ■  0  for  x  <  0  , 

2.  R(x)  exists, 

3.  R(x)  exists, 

4.  R(x)  is  continuous,  and 

5.  R(x)  is  continuous. 

To  prove  (l),  we  note  that  if  R(x)  *  0  for  x  <  0  ,  then  K(x,y)  « 
0  for  x  <  -y  and,  consequently,  V(x)  *  0  for  x  *  0  .  The  proof  of  the 
remaining  four  conditions  requires  somewhat  more  detail. 

Our  starting  point  will  be  the  integral  equation 


0  *  R(x+y)  +  K(x,y)  +  J  R(y+z)K(x,z)dz 


(A-60) 


If  R(x)  is  continuous,  then  Eq.  (A-60)  can  be  written 


0  »  R(x+y)  +  K(x,y)  +  /  R(y+z  )K(x,  z  )dz  .  (A-61) 

-y 


If  R(x)  exists,  Eq.  (A-61)  can  be  differentiated  with  respect  to  x  to 
obtain 


x 

0  *  R(x+y)  +  I^(x,y)  +  K(x,x)R(x+y)  +  f  KK(x/t)R(y+t)dz  .  (A-62) 

-y 


If  R(x)  exists,  then  Eq.  (A-62)  can  be  differentiated  with  respect  to  x  to 
obtain 
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0  -  B(x+y)  +  K^x^y)  +  K(x,x)R(x+y)  +  R(x+y)  ^  K(x,x) 


A 

Kx(x,x)R(x+y)  +  f  K^foyjRfy+zJdz  .  (A-63) 


Next,  differentiate  Eq.  (A-61)  with  respect  to  y  to  obtain 


x 

0  *  R(x+y)  +  K(x,y)  +  K(x,-y)R(0)  +  f  K(x,z)R(y+z)dz 


(A -64 ) 


If  R(x)  is  continuous  and  if  R(x)  *  0  for  x  <  0  ,  then  R(0)  «  0  and 
Eq.  (A-64 )  becomes 


0  *  R(x+y)  +  Ky(x,y)  +  P  K(x,z)R(y+z)dz  .  (A-65) 


Next,  differentiate  Eq.  (A-65)  with  respect  to  y  to  obtain 


A 

0  *  R(x+y)  +  K^y(x,y)  +  K(x,-y)R(0)  +  J  K(x,z)R(y+z)dz  .  (A-66) 


If  R(x)  is  continuous,  then  R(0)  ■  0  and  Eq.  (A-66)  becomes 


A 

0  =  R(x+y)  +  Kyy(x,y)  +  J  K(x,  z  )R(y+z  )dz  .  (A-67) 


Next,  the  integral  in  Eq.  (A-67)  is  integrated  twice  by  parts  and  Eq.  (A-67) 
becomes 


0  »  R(x+y)  +  Kyy(x,y)  +  K(x,x)fl(x+y)  -  Ky(x,x)R(x+y) 

x 

♦/  Kjry(x^*)R(3r+*)4*  • 


(A-66) 
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Next,  Eq.  (A-68)  is  subtracted  from  Eq.  (A-63)  and,  from  this  result,  the 
product  of  2  -  K(x,x)  and  Eq.  (A-61)  is  subtracted.  This  yields 


0  a  "  Kyyfo*)  "  2K(x'y*  ^  K<x>x) 


f  R(y+a)[lC)0C(x,z)-Kyy(x,z)-2K(x,z) K(x,x)]dz  .  (A-69) 

-y 


Obviously,  Eq.  (A-69)  is  satisfied  if 

K^^y)  -  Kyy(x,y)  -  2K(x,y)  ^  K(x,x)  -  0  .  (A-70) 


Subtracting  Eq.  (A-70)  frco  Eq.  (A-54),  we  obtain  the  desired  result: 


V(x)  »2  j-  K(x,x) 

OX 


(A-71) 


Thus,  the  sufficiency  of  conditions  (l)  to  (5)  has  been  demonstrated  and  we 
see  that  if  r(k)  is  a  function  whose  transform  B(x)  satisfies  these  con¬ 
ditions  (as  well  as  the  conditions  specified  in  II),  then  the  construction 
will  be  successful. 


IV.  The  Determination  of  n(x) 

Now  that  we  have  determined  the  scattering  potential,  V(z)  ,  which 
appears  in  the  equation 


u(k,z)  +  [l^-VU)^  u(k,z)  *  0  ,  (A-72 ) 

our  next  task  will  be  to  find  the  index - of -re fr act i on  profile,  n(x)  ,  which 
appears  in  the  equation 


u(k,x)  +  k2n2(x)w(k,x)  »  0 


(A-73) 


The  requirements  imposed  on  V(z)  , 


V(z )  =  0  for  z  <  0 

and 

(A-74) 

lim  V(z)  *  0  , 

z  >°° 

and  those  which  apply  to  n(x)  , 

(A-75 ) 

n(x)  *  1  for  x  <  0 

and 

(A-76) 

lim  n(x)  =  1  , 

x— 

are  sufficient  to  insure  that 

(A-77 ) 

u(k,x)  »  w(k,x)  for  x  <  0 

find 

(A-78) 

lim  j\i(k,x)-w(k,x)]  *  0 

(A-79) 

x— 


Thus,  EqB.  (A-78)  and  (A-79 )  serve  to  prescribe  boundary  conditions  for  Eq. 
(A-73  ). 

In  Eq.  (A-72),  let  x  =  x(z )  .  Then  we  obtain 

u(k,x)  +  u(k,x)  +  — u(k,x)  -  0  (A-80) 

(if  (if 


where  the  dots  indicate  differentiation  with  respect  to  x  and  the  primes 
indicate  differentiation  with  respect  to  z  .  Next,  in  Eq.  (A-80),  let 
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(A-81) 


Then  we  obtain 


u(k,x)  -  v(k,x)(x)-^  . 


v(k,x)  + 


_*L  +  likt  .  JL  .  Viiil  v 

L(x)2  ^(x)2  2(x)3  (xfm 


(k>x)  -  0 


(A-82) 


By  inspection,  ve  see  that  Eq.  (A-83)  is  equivalent  to  Eq.  (A-73)  if 


I 

X 


n(x) 


(A-83) 


and 


(if-® 


-  4V(x)  “  0 


(A-Q4  ) 


Consequently,  if  we  can  find  an  expression  for  x  from  Eq.  (A-84),  we  can 
determine  n(x)  from  Eq.  (A-83).  It  is  also  significant  to  note  that,  because 
n(x)  =  1  for  x  <  1  ,  v(k,x)  in  Eq.  (A-82)  satisfies  the  prescribed  boundary 
conditions . 


Then 


In  Eq.  (A-84),  let 


M 


(A-85) 


(A-86) 


and  Eq.  (A-84)  becomes 


2^(z)  -  y2(z)  +  4V(z)  -  0  .  (A-87) 
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Equation  (A-87 )  will  be  recognized  as  the  Ricatti  equation.  The  substitution 


y(z)  - 


(A-88) 


transforms  it  into 


U(z)  -  V(z)U(z)  «  0 


(A-89 ) 


If  we  solve  Eq.  (A-89)  for  U(z)  ,  we  can  then  find  x  from 


and  x(z)  from 


^(z) 


(A-90) 


(A-91) 


Then,  inverting  x(z)  to  find  z(x)  ,  ve  obtain 


n(x) 


(A-92  ) 


For  x  <  0  ,  n(x)  ■  1  .  Thus  z(x)  *  x  and  U(z)  *  1  for  i  <  0  . 
By  requiring  that  U(z)  and  6(z)  be  continuous,  ve  obtain 

U(0)  -  1  and  U(0)  -  0  .  (A-95) 


These  are  the  boundary  conditions  which  allow  us  to  solve  Eq.  (A-89 )  and  this 
completes  the  problem. 

In  summary,  ve  are  given  the  function  r(k)  and  want  to  find  n(x)  . 
To  do  this,  ve  oust  complete  the  following  sequence  of  operations: 
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1.  Compute  R(x)  from 


r(x)  -  i-  f  r(k)e“ikx  dk 
2tt  u 


(A-94 ) 


where  c  is  a  contour  extending  from  k  ■  -«  to  k  »  +»  and  lying  in  the 
upper  half  of  the  K-plane  above  the  poles  of  r(k)  . 

2.  Find  K(x,y)  from  the  integral  equation 


A 

0  a  R(x+y)  +  K(x,y)  +  J  R(y+z )K(x,z )dz 


(A-95 ) 


3.  Find  V(x)  from 

V(x)  *  2  i  K(x,x)  .  (A-96) 

dx 

4.  Find  U(z)  from 

U(z)  -  V(z)U(z)  -  0  ; 

U(O)  -  1  and  U(o)  «  0  . 

5.  Find  x(z)  from 

x(z) '  ‘  (A'M) 

6.  Invert  x(z)  to  obtain  z(x)  . 

7.  Compute  n(x)  from 

n(x)  -  ^(x)  .  (A-lOO) 

The  constant  c  in  Eq.  (A-99)  can  be  evaluated  if  one  value  of  n(x)  ,  i.e., 
n(0)  la  known. 


(A-97) 
(A-98 ) 
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APPENDIX  B 


THE  MEASUREMENT  OF  THE  REFLECTION  COEFFICIENT.  r(k) 

The  mathematical  theory  for  the  remote  measurement  of  the  index  of 
refraction,  which  was  advanced  earlier  in  this  report,  depended  succinctly 
upon  a  certain  quantity,  r(k)  ,  called  the  reflection  coefficient.  In  this 
section,  ve  consider  the  manner  in  which  this  quantity  may  he  measured. 

The  reflection  coefficient,  r(k)  ,  first  appeared  in  the  equation 

u(k,x)  *  eikx  +  r(k)e"itot  ,  (B-l) 

where  u(k,x)  was  assumed  to  be  a  solution  of  the  (separated)  wave  equation 

u(k,x)  +  k?u(k,x)  a  0  .  (B-2) 

Equation  (B-2)  had,  in  turn,  been  obtained  frcm  the  (total)  wave  equation 

2  2 

u(k,x,t)  -  —  u(k,x,t)  *  0  (B-S) 

dx2  c2  dt2 

by  separation  of  the  temporal  factor  exp(-lkct)  j  that  is,  by  assuming 
that 

u(k,x,t)  ■  u(k,x)e"lkct  .  (B-4) 

Thus,  we  can  write 

u(k,x,t)  -  eilc^x“ct^  +  r(k)e”ilt^x+ct^  .  (B-5) 
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Nov,  u(k,x,t)  as  defined  by  Eq.  (B-5)  represents  a  particular  solution  of 
Eq.  (B-3 ) ;  a  general  solution  can  be  obtained  by  multiplying  by  f(k)  ,  an 
arbitrary  function  of  the  wave  number  k  only.  Let 


g(k)  *  r(k)f (k) 


(B-6) 


and  this  general  solution  is 


uff(k,x,t)  -  f(k)elk(x"ct)  +  g(k)e'ik(x+ct)  .  (B-7) 

6 


This  general  solution  represents  the  sum  of  tvo  waves,  a  forward  (toward 
x  =  +®  )  traveling  wave,  f(k)  exp[ik(x-ct)]  ,  and  a  backward  (toward  x  *  -<*>  ) 
traveling  wave,  g(k)  exp[-ik(x+ct)J  .  Both  waves  travel  with  a  velocity  c 


c 


1 

(u0e0)^ 


(B-8) 


By  further  inspection,  we  see  that  the  forward  traveling  wave  is 
actually  the  sum  of  cm  infinite  number  of  sinusoidal  components,  each  component 
having  a  spatial  frequency,  or  wave  number,  k  ,  and  an  amplitude  of  f  (k)  .  A 
similar  interpretation  can  be  given  to  the  backward  traveling  wave.  Accord¬ 
ingly,  the  total  waves  can  be  obtained  by  summation: 


F(x-ct)  -  -7—  f  f(k)eik(*‘ct)  dk 

V2tt  4m 

G(x+ct)  -  — —  f  g(k)e"ik^x'ct^  dk 


(B-9) 


(B-10) 


Now,  the  two  functions  F(x-ct)  and  G(x+ct )  are  the  space-time 
pr^lles  of  the  incident  and  reflected  waves,  respectively.  They  represent 
two  measurable  quantities,  and  it  is  these  two  quantities  which  we  shall 
measure. 
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Nov  notice  that  F(x-ct)  and  G(x+ct )  are  Fourier  transforms  of 
the  functions  f(k)  and  g(k)  ,  respectively.  Accordingly,  once  F(x-ct) 
and  G(x+ct)  have  been  measured,  f(k)  and  g(k)  can  be  obtained  by  in¬ 
version: 


CO 

f(k)  *  f  F(x-ct)e"ik(x'ct)  d(x-ct)  (B-ll) 

‘ L 

and 

00 

g(k)  *  — —  f  G(x+ct  )e*k(x+at )  d(x+ct)  .  (B-12) 

V2n  L 


Then,  r(k)  may  be  found  from 


r(k) 


f(lc) 


(B-13) 


In  practice,  a  very  short  impulse,  F(x-ct)  ,  will  be  transmitted 
into  the  region  of  nonconstant  index  of  refraction  and  a,  presumably  longer, 
impulse  G(x+ct)  will  be  reflected.  The  requirement  that  F(x-ct)  be  a  very 
short  impulse  stems  from  the  fact  that  the  index-of -refraction  profile  will  be 
slowly  changing  with  time.  If  F(x-ct)  is  too  long,  or,  in  other  words,  if 
too  much  time  is  taken  to  make  the  measurement,  then  a  "smearing”  will  take 
place  which  will  invalidate  the  results,  Rrcm  another  point  of  view,  we  wish 
to  measure  the  value  of  the  reflection  coefficient  over  as  broad  a  band  of 
spatial  frequencies,  k  ,  as  is  practicable.  To  do  this,  it  is  necessary  that 
the  s  pec  trim  of  F(x-ct)  be  as  broad  as  possible  and  one  way  of  obtaining  a 
broad  spectrum  is  to  use  a  very  short  Impulse. 
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APPENDIX  C 


AN  ALTERNATE  DERIVATION  OF  THE  INVERSE  SCATTERING  THECBY 


In  this  section  ve  shall  present  an  alternate  derivation  of  the 
inverse  scattering  theory  and  obtain,  as  a  consequence,  a  different  (but 
equivalent)  inversion  formula  which  may  be  useful  under  sane  circumstances. 
We  shall  also  present  examples  to  demonstrate  the  utility  of  both  inversion 
formulae. 


Formally,  ve  are  given  the  differential  equation 


u(k,x)  +  [k^-Vfc)]  u(k,x)  *  0 


(C-l) 


where 


V(x)  *  0  for  x  <  0 


(C-2) 


and  u(k,x)  is  known  for  x  s  0  .  We  wish  to  find  the  value  of  V(x)  for 
x  2  0  .  New,  assume  that 


u(k,x)  «  f(k)eikx  +  g(k)e"ikx  ,  x  s  0  .  (C-3) 

Following  the  same  procedure  as  in  Appendix  A,  ve  define  the  two  unperturbed 
solutions 


v(k, x)  -  ffrje1**  +  g(k)e-ita  (C-4) 

and 

v(k,x)  *  f#(k)e"ilac  -  g^k)*1**  .  (C-5) 
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Next,  we  define  the  function  |h(k)|2  by 


|h(k)|2  =  |f(k)|2  -  |G(k)|2  .  (C-6 ) 

Consequently,  v(k,x)  has  been  normalized  so  that 


New,  assume  that  there  exists  a  function  K(x,y)  such  that 

L.  r  ^  -  6(x.y)  +  K(x,y)  .  (c-e) 

2n  J-~  lh(k)l2 

The  existence  of  this  function,  K(x,y)  ,  may  be  verified  by  exactly  the  same 
method  used  in  Appendix  A,  Eqs.  (A-47)  to  (A-56),  and  we  find  that  it  must 
satisfy  the  equations 


K^foy)  -  Kyy(x,y )  -  V(x)K(x,y)  «=  0  , 

(C-9) 

£  K(x,x)  =  £  V(x) 

(c-io) 

K(x,-x)  ■  constant 

(c-u) 

By  the  same  token,  using  the  methods  of  Appendix  A,  Eqs.  (A -42)  to  (A-46),  we 
see  that 


Now,  apply  the  operator 


(C-13 ) 


to  Eq.  (C-12)  to  obtain 

x  x 

0  ■  G(y+x)  +  J  K(x, z  )F(y-z )dz  +  J  K(x,z)G(y+z)dz  ,  (C-14) 


-x 


where 


and 


F(x)  *  i-  f  f (k)e"ikx  dk 
2tt  { 


(C-15 ) 


G(x)  mbf  g(k)e’ikx  ** 

c 


(C-16 ) 


Equation  { C-14 )  is  the  alternate  Inversion  formula. 

To  indicate  the  equivalence  of  the  two  inversion  formulae,  as  well 
as  to  demonstrate  one  method  of  solution,  ve  shall  consider  two  examples .  In 
the  first  example,  let 


and 


f(k) 


l(k-i) 

(k+i)2 


g(k) 


(k+i)3 


(C-17 ) 


(C-18) 


Then 


r(k) 


g(k)  „  -1 

f(k)  n2+1 


(C-19) 
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Using  the  calculus  of  residues, ve  may  easily  ascertain  that 


F(x)  «  (l-2x)e‘x  ,  x  i  0 


F(x)  “0  ,  x  <  0  , 

(C-20) 

G(x)  *  ^  x2e”x  ,  x  2  0 

G(x)  ■  0  ,  x  s  0  , 

(C-21) 

and  that 

R(x)  »  sinh  x  ,  x  i  0 

R(x)  =  0  ,  x  s  0 

(C-22) 

Consequently,  we  are 

concerned  with  the  two  integral  equations 

0  =  G(x+y) 

+  J  K(x,z )F(y-z  )dz  +  /*  K(x,  z  )G(y+z  )dz 
-x  -y 

(C-23) 

and 

0  ■  R(x+y)  +  K(x,y)  +  J  K(x,z  )R(y+z  )dz 

-y 

(C-24) 

We  shall  solve  Eq.  (C-24)  first. 

Note  that 

HvMy)  •  o 

(C-2S) 

where 

p(p)  *  p2-! 

(C-26) 
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and 


P 


Applying  this  operator  to  Eq.  (C-24),  we  ctotain 


0  -  P(p)R(x+y)  +  P(p)K(x,y)  +  f  K(x,z )P(p)R(y+z )dz 

-y 


+  K(x,-y)  [r(0)+R(0)]  -  Ky(x,-y)R(0)  . 

Now,  because  R(0)  =  0  ,  and  R(0)  =  1  ,  Eq.  (C-28)  becomes 

0  -  K^Ujy)  -  K(x,y)  +  K(x,-y)  . 

The  solution  of  Eq.  (C-29)  is 

K(x,y)  -  a(x)  +  b(x)  sinhViy 


(C-27) 


(C-28) 


(C-29) 


(C-30) 


where  a(x)  and  b(x)  are  arbitrary  functions  of  x  alone.  By  substituting 
the  expression  for  K(x,y)  from  Eq.  (C-30)  into  Eq.  (C-24 )  and  preforming  the 
indicated  integrations,  we  obtain 


K(x,y) 


1  sinh  Vix  +  sinh  Viy 
— =  . '.r  . ,  y  i  ~x 

Vi  cosh  V5x 


K(x,y)  -  0 


,  y  <  -x 


(C-31) 


Consequently,  from  Eq.(C-lO)  we  easily  obtain 

V(x)  -  -4  sinh2  Vix  ,  x  i  0 

V(x)  ■  0  ,  x  <  0  .  (C-32) 
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Next,  ve  turn  cur  attention  to  Eq.  (C-23).  We  obtain  easily 


P(p)F(y)  -  0 


(C-33) 


where 


P(p)  »»  +  2p  +  1 


(C-34) 


and  that 


(p+l)P(p)G(y)  -  0  .  (C-35) 

Thus,  by  applying  the  operator  (p+l)3  to  Eq.  (C-23),  and  using  the  values 

A(0)  -  1  ,  A(0)  *  -3  ,  A(0)  -  5  (C-36 ) 

and 

B(0)  -  0  ,  B(0)  “  0  ,  B(0)  »  1  ,  (C-37 ) 


we  obtain 


0  -  Kyy(x,y)  -  K(x,y)  +  K(x,-y)  .  (C-38) 


Equation  (C-38)  is  identical  to  Eq.  (C-29).  Consequently,  its  solution  is 
given  by  Eq.  (C-31)  and  V(x)  by  Eq.  (C-32). 

It  is  important  to  note  frcm  this  example  that,  while  R(x)  and 
G(x)  are  continuous,  A(x)  has  a  discontinuity  at  the  origin.  Frcn  this,  we 
may  assume  that  it  is  not  necessary  that  A(x)  be  continuous.  However,  it 
would  be  incorrect  to  assume,  in  this  respect,  that  Eq.  (C-23)  represents  a 
generalization  of  Eq.  (C-24 )  as  the  next  example  will  show. 
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For  the  second  example,  let 


f(k) 


l 

k+i 


sOO 


l 

(k+i )2 


and 


r 


00  - 

f(k) 


-i 

k+i 


Then 


R(x)  =  -e'x  ,  x  a  0 
=  0  ,  x  <  0 


(C-39) 


(C-40) 


(C-41) 


In  this  case,  R(x)  is  discontinuous  and  we  would  expect  the  theory  to  fail 
(which  it  does).  R(x)  satisfies 


(p+l)R(x)  =«  0  .  (C-42) 


Thus,  by  inserting  the  value  of  R(x)  into  Eq.  (C-24 )  and  performing  the 
operation  (p+l)  ,  we  obtain 


0  ■  Ky(x,y)  +  K(x,y)  -  K(x,-y)  .  (C-43) 


The  solution  of  Eq.  (C-43)  is 


K(x,y)  •  K(x) 


(C-44) 


But,  from  the  relation 


K(x, -x)  ■  constant 


(C-45) 


•  58  • 


ve  obtain 


K(x,x)  »  constant  (c-46) 

and 

V(x)  =  0  ,  all  x  .  (C-47) 

Hcvever,  the  reflection  coefficient  r(k)  *  0  corresponds  to  V(x)  *  0  ,  not 
the  one  given  in  Eq.  (C-40). 

Continuing  the  example,  ve  have 

F(x)  *  e*x  ,  x  a  0 

F(x)  =0  ,  x  <  0  (C-48) 

and 

G(x)  *  -xe"x  ,  x  i  0 

G(x)  “0  ,  x  £  0  .  (C-49) 


These  functions  satisfy  the  equation 


(p+l)F(x)  -  (p+l )^G(x)  -  0 


(C-50) 


Thus,  by  substituting  the  value  of  F(x)  and  G(x)  into  Eq.  (C-23)  and  per¬ 
forming  the  operation  (p+1 ve  obtain 


0  -  Ky(x,y)  +  K(x,y)  -  K(x,-y) 


(C-Sl) 


and,  as  before,  the  theory  fails. 

In  concluding  this  Appendix,  let  us  point  out  that  the  second  exam¬ 
ple,  used  here  to  demonstrate  the  failure  of  the  theory,  is  exactly  the  same 
as  the  example  used  In  Appendix  C  of  the  Fourteenth  Interim  Technical  Report 
(10  July  1962)  to  demonstrate  the  validity  of  the  theory.  The  answer  to  this 
unseemly  contradiction  is  that  all  the  work  appearing  in  Appendix  C  of  the 
Fourteenth  Interim  Technical  Report  is  incorrect. 
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APPENDIX  D 


CALIBRATED  TRANSMISSION  PATH 


The  calibrated  transmission  path  will  consist  of  a  section  of  K 
band  waveguide  (either  RG-91/U  brass  guide  or  RG-107/U  silver  guide)  approxi¬ 
mately  one  meter  long.  This  waveguide  will  be  carefully  assembled  and  cali¬ 
brated  and  will  provide  the  experimenter  with  a  transmission  path  whose  index- 
of-refraction  profile  and  other  pertinent  characteristics  are  accurately 
known.  Initially,  experiments  will  be  performed  to  indicate  how  variations  in 
the  index-of -refraction  profile  along  the  line  will  affect  length  measurements 
made  using  a  phase-null  or  interferometer  technique.  Later  experiments  will 
be  aimed  at  determining  how  accurately  the  reflections  caused  by  the  varia¬ 
tions  of  the  index-of -refraction  profile  can  be  measured  and  how  these  meas¬ 
urements  can  be  used  to  more  accurately  determine  the  length  of  the  path. 

The  length  of  the  guide  will  be  determined  using  the  equipment  shown 
in  Fig.  6.  This  equipment  will  consist  of  the  following  components: 

1.  Huggins  Model  615  Oscillator  (10-20  Gc.,  0.5  mw.), 

2.  Huggins  Model  322  TWT  Amplifier  (30  db,  1  w.), 

3.  Waveguide  Tee  ( Demo may  Bonardi  DBF620), 

4.  Crossguide  Coupler  (DBF 631 ), 

5.  Matched  Load  (DBF450), 

6.  Elbow  (DBF224), 

7.  Attenuator,  0-50  db  (DBF420), 

8.  Thermistor  Mount  (Hewlett-Packard  F467B),  and 

9.  Power  Meter  (Hewlett-Packard  430  C). 

To  determine  the  waveguide  length  using  the  phase-null  system,  a 
short  circuit  will  be  placed  at  one  end  of  the  guide  and  a  microwave  signal 
from  the  phase-null  equipment,  Fig.  6,  fed  into  the  other.  The  electric  field 
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within  the  guide  will  be  zero  at  the  short  circuited  end  and,  because  of  the 
standing  waves  caused  by  the  short,  will  also  be  zero  every  one-half  guide 
wavelength  along  the  line  measured  from  the  short. 

The  oscillator  frequency  will  be  varied  until  the  voltage  at  the 
sending  end  is  zero,  and  the  half  guide  wavelength,  X^  ,  corresponding  to 
this  frequency  will  be  measured.  Then,  the  length  of  the  guide,  L  ,  can  be 
expressed  as 


L  *  mX^ 


where  m  is  an  integer  yet  to  be  determined.  Next,  increase  the  oscillator 
frequency.  The  voltage  at  the  sending  end  will  increase  (as  the  frequency 
increases),  reach  a  maximum,  and  then  decrease  to  zero.  Denote  the  half  guide 
wavelength  at  which  this  zero  occurs  as  X2  .  The  waveguide  length  can  now  be 
expressed  in  terms  of  X2  as 

L  = 


From  these  two  results,  it  follows 


L 


and,  since  both  X^  and  Xg  sure  known,  L  can  be  determined. 

Unfortunately,  computational  difficulties  arise  at  this  point  be¬ 
cause  of  the  necessity  of  dividing  by  the  difference  of  two  numbers,  X^  and 
X2  ,  which  are  nearly  equal.  To  illustrate,  assume  that  X^  *  2  cm.  and  that 
L  =  100  cm.  Then  m  ■  50  and 


(m+l)  X2 


that 


XlXg 

xx-x2 


\  .  i22  »  1.960784  ...  cm. 
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Nov,  assume  that  X  can  only  be  measured  to  an  accuracy  of  t 0.001  cm.  Then 
the  two  measured  quantities  would  be 
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X^  *  2.000  cm., 
Xg  =  1.961  cm. 


and 

h  ■  (2-Q00)(1.961)  ,  na.se.  ... 
2t000-9«161 


resulting  in  an  error  of  0.564  cm.  in  L  . 

To  remedy  this  problem,  assume  that,  as  before,  X^  represents  the 
half  guide  wavelength  at  which  a  zero  of  the  input  voltage  occurs.  Then  in¬ 
crease  the  oscillator  frequency  noting  that  the  input  voltage  varies  through 
successive  maxima  and  zeros.  Let  Xn  be  the  half  guide  wavelength  at  which 
the  nth  zero  occurs.  Then 


1  =  mX^  »  (m+n)Xn 


and 


1  =  n 


XlXn 

Xl"Xn 


Now,  assume  that  X^  *  2  cm.,  L  «  100  cm.,  and  n  =  40  .  Then 
X  a  1 . Ill  ...  cm . 

Assume,  once  again,  that  Xn  can  be  measured  to  an  accuracy  of  tO.OOl  cm. 
Then  the  measured  quantities  will  be 


Xx  =  2.000 

X40  “  1-m 
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and 


1  =  40  (g’000)(l«lH)  x  99.978 

2.000-1.111 


or,  an  error  of  only  0.022  cm.  in  L  . 

These  numerical  examples  serve  only  to  Illustrate  the  nature  of  the 
error  which  will  result  and  also  to  Indicate  that  the  error  can  be  minimized 
by  maximizing  the  term  "n"  in  the  formula 


L  = 


n 


XlXn 

Xl‘Xn 


For  a  line  one  meter  long  and  a  range  of  frequencies  12-18  Gc.,the  maximum 
value  of  n  will  be  approximately  50  (when  X^  *  21^  *  2  cm.).  TO  obtain 
an  approximation  of  the  error  which  can  be  expected,  use  the  above  values  to 
write 


AL  »  50 


(Sf*l)(l+e5o) 

(2+ejMl+eso) 


100 


where  AL  is  the  error  in  L  and  and  «5Q  are  the  errors  in  X^^  and 
X^  .  Expanding  the  above  equations. 


AL  =>  50 


tl*50~cl 

l+«r«50 


Next,  assume  that  l«ll  “  I *50  f  =  0.001  cm.  and  maximize  the  numerator  and 
minimize  the  denominator  to  obtain 


AL  <  50 


0.001001 

0.998 


■  0.05  cm. 


In  other  words,  with  the  equipment  specified,  the  length  of  a  one-meter 
section  of  guide  can  be  determined  to  better  than  !i  mm. 
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Before  making  measurements,  the  test  equipment,  Fig.  6,  must  be  cali¬ 
brated.  To  do  this,  connect  a  calibrated -movable- short  such  as  Demomay 
Bonardi  DBF969-1  to  the  test  equipment  (the  open  arm  of  the  tee).  Then, 
follow  the  procedure  below: 

1.  Set  the  oscillator  to  that  frequency  near  12  Gc.  where  zero 
power  is  indicated  (on  the  430  C  power  meter). 

2.  Move  the  short  in  a  direction  away  from  the  tee,  noting  that  the 
power  indicated  on  the  430  C  will  increase  and  then  decrease  again  to  zero. 
Record  the  distance  which  the  short  has  moved  as  X1  and  return  the  short  to 
its  original  position. 

3.  Increase  the  oscillator  frequency  to  about  18  Gc. noting  that  as 
the  frequency  is  increased,  the  indicated  power  varies  through  successive 
maxima  and  zeros .  Set  the  oscillator  to  that  frequency  where  a  zero  is  ob¬ 
tained  and  record  the  total  number  of  zeros  as  n  . 

4.  Measure  this  new  wavelength  by  the  procedure  in  (2)  and  record 
this  value  as  Xn  . 

5.  Compute  the  effective  length  of  the  calibration  equipment,  Le  , 

from 


L 

e 


xxx 

X1-X] 


This  length  must  be  subtracted  from  all  subsequent  measurements.  It  should, 
of  course,  be  measured  several  times  and  the  results  averaged  to  obtain  maxi- 
mum  accuracy. 

To  measure  the  length  of  a  section  of  waveguide,  connect  the  cali¬ 
brated  length  measuring  equipment  to  one  end  of  the  guide  and  the  short  to  the 
other.  Repeat  the  procedure  above  to  determine  the  combined  length  of  the 
guide  and  effective  length  of  the  equipment,  L+Le  .  Subtract  Le  from  the 
result  and  the  remainder  will  be  the  length  of  the  guide. 

Once  the  measuring  equipment  and  test  line  have  been  set  up  and 
calibrated,  the  next  step  will  be  to  Introduce  artificial  dielectric  perturba¬ 
tions  into  the  test  line  and  the  experimenter  will  attempt  to  (a)  measure  the 
reflections  produced  by  the  perturbations,  and  (b)  use  these  measurements  to 
prove  the  validity  of  the  inverse  scattering  theory. 
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Whenever  a  dielectric  is  introduced  into  a  waveguide,  the  effect  will 
he  to  increase  the  length  of  the  guide.  This  fact  can  be  demonstrated  with  a 
simple  example. 

Assume  that  the  guide  is  100  cm.  long  and  that  the  operating  fre¬ 
quency  is  15  Gc.  Then,  for  an  air  filled  guide  (*r  ■  l),  the  guide  wavelength 
will  be  about  2.5  cm.  (refer  to  Pig.  7  which  shows  the  relationship  between 
guide  wavelength  and  frequency  for  guides  filled  with  various  dielectrics). 
Thus,  the  length  of  the  guide  will  be  100  +  2.5  ■  40  wavelengths.  Now, 
suppose  that  a  block  of  dielectric  having  a  relative  dielectric  constant  of 
cr  ■  20  and  a  length  of  10  cm.  is  inserted  into  the  guide.  Referring  to 
Fig.  7,  one  sees  that  the  wavelength  in  this  dielectric  is  about  0.45  cm. 

Thus,  the  line  length  will  be  90  +  2.5  +  10  +  0.45  =  145  wavelengths,  an  in¬ 
crease  of  over  100  wavelengths.  If  one  were  to  attempt  to  measure  the  length 
of  this  line  using  the  phase-null  method  just  described,  the  result  would  be 
145  x  2.5  =  362.5  cm.  It  is  toward  the  correction  of  errors  of  this  nature 
that  the  present  research  is  directed. 

Initially,  the  reflection  coefficient  of  the  artificial  dielectric 
perturbations  will  be  measured  using  standard  techniques  and  the  data  so  ob¬ 
tained  will  be  used  to  prove  the  validity  (and  utility)  of  the  inverse  scat¬ 
tering  theory.  However,  it  is  apparent  that,  in  order  to  measure  reflection 
coefficients  of  -120  db  or  less,  new  techniques  will  have  to  be  perfected. 
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APPENDIX  E 


RADIO  FREQUENCY  ENERGY  SOURCE 


For  actual  index-of -refraction  profile  measurements,  the  tunable, 

C-W  source  used  in  the  calibration  of  the  transmission  path  will  be  replaced 
by  a  pulsed  source.  This  pulsed  source  will  consist  of  a  pulse -producing 
system  and  a  pulse-shaping  system.  The  pulse-producing  system  will  comprise 
an  ordinary  Ku-band  magnetron  and  associated  modulation  circuitry  for  pro¬ 
ducing  pulses  of  rf  energy  (either  on  a  single-shot  or  a  repetitive  basis) 
which  are  of  the  order  of  one  microsecond  in  duration.  The  pulse-shaping  sys¬ 
tem,  at  the  output  of  the  magnetron,  will  convert  these  microsecond-pulses  into 
the  extremely  short-duration  pulses  (of  the  order  of  a  few  nanoseconds)  that 
will  be  necessary  for  the  experiments.  Since  the  development  of  the  pulse¬ 
shaping  system  represents  a  major  area  of  research  effort,  very  little  can  be 
said  about  the  final  form  of  this  portion  of  the  pulsed  source.  However,  the 
design  of  the  pulse-producing  portion  of  the  pulsed  source  (i.e.,  the  modulator 
and  magnetron)  can  be  specified  within  the  existing  state  of  the  art. 

Experience  has  shown  that  a  QK319,  Ku-band  magnetron  can  be  success¬ 
fully  pulsed  with  an  APS-2,  hard- tube  radar  modulator.  Since  both  the  QK319 
magnetron  and  the  APS-2  modulator  are  readily  available,  at  very  reasonable 
cost,  a  pulse-producing  system  based  on  these  items  is  contemplated.  Some 
construction  work  will  be  necessary  in  the  preparation  of  a  modulator  control 
system,  and  an  appropriate  means  for  monitoring  the  operation  of  the  pulse- 
producing  system  will  have  to  be  provided. 

A  block  diagram,  showing  the  essential  features  of  the  pulse-producing 
system,  is  presented  in  Fig.  6.  Typical  features  of  the  APS-2  radar  modulator 
are  shown  in  the  schematic  diagram  of  Fig.  9.  A  triggering  signal,  provided 
through  contact  F  on  J805  is  impressed  on  the  pulse  forming  network,  which  Is 
shown  in  the  lower  right-hand  corner  of  the  diagram.  The  actual  arrangement  of 
the  pulse  forming  network  is  determined  by  two  relays,  K803  and  K804,  which  are 
controlled,  through  contacts  A,  E,  and  H  of  J605,  from  the  "pulse  duration" 
selector  in  the  modulator  control  system.  In  turn,  the  arrangement  of  the 
pulse  forming  network  determines  the  duration  of  the  pulses  produced  by  the 
driver  tube,  V802,  and  consequently  the  duration  of  the  pulses  developed  by  the 
output  tube,  V805.  Pulse  duration  of  0.5,  1.0,  and  2.0  ^sec.  is  available 
in  all  models  of  the  APS-2  modulator. 

Normally,  the  output  pulses  developed  by  V805  are  delivered  directly 
to  the  cathode  of  a  2722  (S-band)  magnetron  of  a  7 25 -A  (x-band)  magnetron  which 
is  mounted  in  the  lower  portion  of  the  APS-2  modulator.  However,  by  means  of 
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Fig.  8  -  Pulse  Producing  System  Block  Mi 


Fi*.  9  -  Typical  Modulator  Circuit  (APS- 2) 


TO 


a  slight  modification  of  the  APS-2  modulator  circuitry,  the  output  pulses  can 
he  delivered  to  the  cathode  of  a  QK319  (Ku-hand)  magnetron  through  a  short 
length  of  high  voltage  cable.  Experience  has  shown  this  arrangement  to  be 
quite  satisfactory  in  that  it  provides  for  a  great  deal  of  flexibility  in  the 
physical  location  of  the  magnetron.  In  other  words,  the  magnetron,  along  with 
its  filament  transformer,  can  be  mounted  at  any  desired  location  in  an  experi¬ 
mental  setup  (or  even  moved  from  me  location  to  another)  without  the  necessity 
of  "udring  cospllcated  changes  in  waveguide  connections  to  the  heavy,  rather 
Immobile,  modulator  unit. 

It  Should  be  noted  that  under  typical  operating  conditions,  the 
QK319  magnetron  requires  a  15-kv. pulse  for  satisfactory  operation.  Since  the 
APS-2  radar  modulator  is  capable  of  producing  only  a  12-kv.  pulse,  a  modifica¬ 
tion  of  the  QK319  magnetic  circuit  is  necessary.  Although  such  a  modification 
results  in  a  reduction  Of  the  power  output  of  the  QK319,  peak  powers  of  as  much 
as  15-20  kw.  can  still  be  obtained  with  this  magnetron-modulator  combination. 

Control  and  triggering  of  the  APS-2  radar  modulator  can  be  accom¬ 
plished  with  the  special  modulator  control  system  shown  schematically  in  Fig. 
10.  This  modulator  control  circuit  includes  a  10  kc,  crystal  controlled  oscil¬ 
lator  which  provides  a  stable  reference  signal  for  pulse  repetition  frequency. 
Five  flip-flop  circuits  are  used  as  frequency  dividers;  a  suitable  switching 
circuit  allows  nine  different  arrangements  of  the  five  flip-flop  circuits  and 
consequently  provides  nine  different  (crystal  controlled)  pulse  repetition 
frequencies  between  approximately  250  pps  and  1,650  pps.  The  output  of  the 
frequency  divider  circuit  drives  a  blocking  oscillator  which  provides  the 
trigger  pulses  for  the  APS-2  radar  modulator  through  contact  F  on  P805.  "On- 
off"  control  of  the  modulator  plate  supply  is  accomplished  with  the  circuitry 
shown  in  the  upper,  right-hand  portion  of  Fig.  10.  The  various  voltages  re¬ 
quired  for  the  modulator  control  system  are  furnished  by  the  power  supply 
shown  schematically  in  Fig.  11. 

A  Tektronix  RML6  oscilloscope  provides  means  for  monitoring  signals 
within  the  pulse  producing  system.  For  convenience  and  flexibility  of  opera¬ 
tion,  the  monitor  oscilloscope,  APS-2  radar  modulator  and  modulator  control 
system  (including  its  power  supply)  can  be  arranged  in  a  portable  modulator 
unit  as  shown  in  Fig.  12. 
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Pig.  11  -  Power  Supply 
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Fig.  12  -  Portable  Radar  Modulator 
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